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ABSTRACT 

The Second Planck Catalogue of Compact Sources is a list of discrete objects detected in single-frequency maps from the full duration of the 
Planck mission and supersedes previous versions. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. 
Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two 
subcatalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these (PCCS2) covers most of the sky and allows the 
user to produce subsamples at higher reliabilities than the target 80 % integral reliability of the catalogue. The second (PCCS2E) contains sources 
detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E 
include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive 
Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects 
in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and the Early Release 
Compact Source Catalogue (ERCSC). 
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1. Introduction 

This paper, one of a set associated with the 2015 release of data 
from the Planck mission,^ describes the second release of the 


^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
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Table 1. Characteristics of the PCCS2 and PCCS2E catalogues: 
flux density at 90% completeness in total intensity; number of 
sources detected in each catalogue in total intensity; and num¬ 
ber of sources that have a polarized signal is measured above a 
99.99% confldence level. See Tables 13 and 14 for more details. 


Channel 

Flux density 90% 
completeness 
[mJy] 

No. of sources 

Polarized sources 

PCCS2 

PCCS2E 

PCCS2 

PCCS2E 

30 ... . 

427 

1560 


122 


44 ... . 

692 

934 


30 


70 ... . 

501 

1296 


34 


100 ... . 

269 

1742 

2487 

20 

43 

143 ... . 

177 

2160 

4139 

25 

111 

217 ... . 

152 

2135 

16842 

11 

325 

353 ... . 

304 

1344 

22665 

1 

666 

545 ... . 

555 

1694 

31068 



857 ... . 

791 

4891 

43290 




Planck Catalogue of Compact Sources (PCCS2)^. It outlines the 
construction of the single-frequency catalogues from an analy¬ 
sis of each of the nine Planck frequency-channel, full-mission 
maps. The construction of these catalogues builds on much of 
the same infrastructure and methodology as the first incarnation 
of the Planck Catalogue of Compact Sources (PCCS), and the 
reader is referred to Planck Collaboration XXVIII (2014) for a 
full description of the catalogue construction procedures, which 
are only summarized here. Table 1 lists some basic properties of 
the nine frequency subcatalogues of the current release. 

One of the primary differences of this release from the PCCS 
is the division of the six highest frequency catalogues into two 
subcatalogues, the PCCS2 and the PCCS2E. This division sep¬ 
arates sources for which the reliability (the fraction of sources 
above a given S/N which are real) can be quantified (PCCS2) 
from those of unknown reliability (PCCS2E). This separation 
is primarily based on the Galactic coordinates of the source, as 
described in Sect. 2.3. The target integral reliability of the the 
entire catalogue, as in the PCCS, is 80 % or greater. The advan¬ 
tage of setting the reliability target this low is that it improves 
the odds of discovering interesting sources with unusual prop¬ 
erties, which might otherwise have been rejected by restrictive 
selection criteria. On the other hand, a highly reliable catalogue 
is desirable for follow-up observations. That is the aim of the 
PCCS2. To this end, we have provided additional information 
in the catalogue which will allow the user to select a subset of 
highly reliable sources from the PCCS2. This takes the form of 
an additional flag per source that indicates the highest reliabil¬ 
ity catalogue to which that source belongs, allowing the user to 
perform a cut on the PCCS2 to reduce it to the desired percent¬ 
age reliability subset. To assist users, we also flag those sources 
identified in other catalogues, mainly at radio wavelengths. 

The principal data-driven difference between the PCCS and 
the PCCS2 and PCCS2E catalogues is the additional data from 
the “extended” mission and the inclusion of polarization mea¬ 
surements in seven out of the nine frequency channels. The 
545 and 857 GHz channels are not sensitive to polarized sig- 


entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 

^ http://archives.esac.esa.int/pla 
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Fig. 1. Sensitivity (the flux density at 90 % completeness) of the 
PCCS2, compared with PCCS, ERCSC, WMAP and others as 
described in text. The sensitivities displayed for the LEI channels 
are for the full sky. Eor the HEI channels, the 90 % completeness 
limits plotted for the PCCS were evaluated in the extragalactic 
zone, as defined in Planck Collaboration XXVIII (2014). The 
regions of sky to which the 90 % completeness limits apply are 
therefore similar (but not identical) to those of the PCCS2. These 
comparisons are discussed further in Sect. 4. 


nals, so the polarization measurements span the range from 30 to 
353 GHz. The polarization measurements provided are based on 
the positions of the compact sources discovered in the temper¬ 
ature maps; there is no independent search for compact sources 
in polarization. The additional data, together with improved data 
processing, have the effect of reducing the noise and hence im¬ 
proving the completeness of the PCCS2 and PCCS2E catalogues 
over that of the PCCS, as demonstrated in Pig. 1. In this fig¬ 
ure we compare the sensitivity of the PCCS2 with that of the 
PCCS (Planck Collaboration XXVIII 2014), the Early Release 
Compact Source Catalogue (see Planck Collaboration XIII2011 
for the 30-70 GHz channels and Planck Collaboration Int. VII 
2013 for the others), WMAP (see Gonzalez-Nuevo et al. 2008 
for the channels at 41 GHz and below, and Lanz et al. 2013 for 
61 and 94 GHz). The sensitivity levels for the Herschel SPIRE 
and PACS instruments are from Clements et al. (2010) and Rigby 
et al. (201 1), respectively. The other wide-area surveys shown as 
a comparison are: the Green Bank 6 cm Survey, GB6 (Gregory 
et al. 1996), the Combined Radio All-Sky Targeted Eight GHz 
Survey, CRATES (Healey et al. 2007), the Australia Telescope 
20 GHz Survey, AT20G (Murphy et al. 2010), the Planck-ATCA 
Co-eval Observation Project, PACO (Bonavera et al. 2011), 
the South Pole Telescope, SPT (Mocanu et al. 2013a), the 
Atacama Cosmology Telescope, ACT (Marsden et al. 2013), and 
the Infrared Astronomical Satellite catalogue, IRAS (Beichman 
et al. 1988). 

In Sect. 2 we describe the construction of the catalogues, in¬ 
cluding the criteria used to place a source into the PCCS2 or 
the PCCS2E, and the methods used to measure the flux densi¬ 
ties and linear polarization parameters. Section 3 discusses the 
validation and quality assessment of the catalogues. Here the 
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internal and external consistency tests are described, as are the 
completeness (the fraction of the total number of sources above 
a given flux density that are present in the catalogue) and relia¬ 
bility of the catalogues. The overall characteristics of the PCCS2 
and PCCS2E are presented in detail in Sect. 4; they are also com¬ 
pared with the characteristics of the PCCS. The released product, 
which is composed of the catalogues and their associated maps, 
is described in Sect. 5. We summarize our conclusions in Sect. 6. 
Details of the estimators used for photometry and for polariza¬ 
tion measurements are given in the appendices. 

2, The Second Planck Catalogue of Compact 
Sources 

2.1. Data 

After four years of operations, the data from the full Planck 
mission have been transformed into full-sky HEALPix^ maps 
(Gorski et al. 2005) by the data processing centres (DPCs) 
(Planck Collaboration VI 2016; Planck Collaboration VIII 
2016). The Low Frequency Instrument (LFI) DPC produced the 
30, 44, and 70 GHz maps after the completion of eight full sur¬ 
veys (spanning the period 12 August 2009 to 3 August 2013). 
In addition, special LFI maps covering the period 1 April 2013 
to 30 June 2013 were produced in order to compare the Planck 
flux-density scales with those of the Very Large Array and the 
Australia Telescope Compact Array, by performing simultane¬ 
ous observations of a sample of sources over that period. The 
High Frequency Instrument (HFI) DPC produced the 100, 143, 
217, 353, 545, and 857 GHz maps after flve full surveys (2009 
August 12 to 2012 January 11). The flux densities of all sources 
measured in the full-mission maps are an average of several ob¬ 
servations. Single-survey maps are available from the Planck 
Legacy Archive,"^ but single-survey flux densities are not pro¬ 
vided in this catalogue. It is important to note that even single¬ 
survey maps may include more than one observation of an in¬ 
dividual source, and hence extracting flux densities from the 
single-survey maps does not guarantee a single-epoch observa¬ 
tion for a given source. However, in the Planck Legacy Archive 
the time-ordered data are available and users can produce maps 
from arbitrary time intervals. Table 2 gives the parameters of the 
Planck beams used in this paper; further details may be found 
in Planck Collaboration II (2016) and Planck Collaboration VII 
(2016). 

2.2. Catalogue construction 

The compact sources in the catalogue were detected at each fre¬ 
quency independently using improved versions of the detection 
pipelines used to create the PCCS. These pipelines are based on 
the Mexican Hat Wavelet 2 algorithm (MHW2; Gonzalez-Nuevo 
et al. 2006; Lopez-Caniego et al. 2006). This is a cleaning and 
denoising algorithm used to convolve the maps, preserving the 
amplitude of the sources while greatly reducing the large-scale 
structures visible at these frequencies (e.g. diffuse Galactic emis¬ 
sion) and small scale fluctuations (e.g. instrumental noise) in the 
vicinity of the sources. The LFI and HFI DPCs have different im¬ 
plementations of the MHW2 algorithm, but consistency checks 
between them have been performed. The differences between the 
implementations are described in Planck Collaboration XXVIII 
(2014), and are due to the different characteristics of the maps 

^ http://healpix.sourceforge.net 
http: //archives. esac. esa. int/pla 


Table 2. Parameters of the beams used in the construction of the 
PCCS2 and PCCS2E. 


Channel 

FWHM 


Beam area 


Fitted 

[arcmin] 

Effective 

[arcmin] 

Dbeam 

[arcmin^] 

Dbeamj 

[arcmin^] 

Dbeam2 

[arcmin^] 

30 ... . 

32.29 

32.41 

1190.06 

1117.30 

1188.93 

44 ... . 

27.00 

27.10 

832.00 

758.00 

832.00 

70 ... . 

13.21 

13.32 

200.90 

186.10 

200.59 

100 ... . 

9.66 

9.69 

106.22 

100.78 

106.03 

143 ... . 

7.22 

7.30 

60.44 

56.97 

60.21 

217 ... . 

4.90 

5.02 

28.57 

26.46 

28.46 

353 ... . 

4.92 

4.94 

27.69 

25.32 

27.53 

545 ... . 

4.68 

4.83 

26.44 

24.06 

26.09 

857 ... . 

4.22 

4.64 

24.37 

22.58 

23.93 


Notes. Two FWHM values are given: one from an elliptical Gaussian 
fit to the beam, and another that is the FWHM of a Gaussian with the 
same solid angle as the main beam, Dbeam- The FWHM found from the 
main beam solid angle is used to evaluate Dbeami and Dbeam 2 ^ which are 
the beam solid angles within a radius equal to this FWHM and twice 
this FWHM, respectively. 

at LFI and HFI frequencies, requiring alternative methods to re¬ 
duce the numbers of spurious detections. Both implementations 
project the full-sky maps onto square patches where the Alter¬ 
ing and detection is performed. The sizes of the patches and the 
overlap between patches have been chosen in such a way that 
the full sky is effectively covered. Sources above a fixed signal- 
to-noise ratio (S/N) threshold are selected and their positions are 
translated from patch to spherical sky coordinates. Because the 
patches overlap, multiple detections of the same object can oc¬ 
cur; these must be found and removed, keeping the detection 
with the highest S/N for inclusion in the catalogue. 

2.3. Defining the PCCS2 and PCCS2E 

For reasons explained in Sect. 1, for this release we provide 
two subcatalogues for some frequency channels, which we call 
PCCS2 and PCCS2E. We also provide a new parameter for each 
source that gives the highest reliability catalogue to which the 
source belongs. However, it is not possible to evaluate this pa¬ 
rameter for every source. It is this consideration which separates 
the sources into the PCCS2 and the PCCS2E; those sources with¬ 
out this held are placed in the PCCS2E. Figure 2 shows the dis¬ 
tribution of the PCCS2 and PCCS2E sources across the sky for 
three of the Planck frequency channels (30, 143, and 857 GHz). 

LFI: One measure of the reliability of the sources detected in 
the three LFI channels is based on a comparison with existing 
catalogues of radio sources. First, each single-frequency cat¬ 
alogue is compared with the appropriate external radio cata¬ 
logues. Next, all the remaining unidentified sources outside a 
\b\ > 20° Galactic cut are examined on a source-by-source ba¬ 
sis by performing a manual search in archival repositories such 
as the NASA/IPAC Extragalactic Database (NED). Second, a 
multifrequency analysis is used to assess whether or not a source 
is present in more than one Planck channel between 30 and 
100 GHz. All sources with plausible identifications in external 
catalogues are assigned to the PCCS2 with a high degree of reli¬ 
ability. Additionally, all sources that are detected in two or more 
Planck channels in the multifrequency analysis are also placed 
in the PCCS2, albeit with a lower degree of reliability. Given the 
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Fig. 2. Top: distribution of the validated sources from the PCCS2. Red, blue, and green circles show sources from the 30, 143, and 
857 GHz catalogues, respectively. Bottom: locations of the sources in the 143 and 857 PCCS2E, shown by blue and green circles, 
respectively. The figure is a full-sky Mollweide projection with the Galactic equator horizontal; longitude increases to the left with 
the Galactic centre in the centre of the map. The size of the filled circles gives an idea of the relative flux densities of the sources 
per frequency, where the larger circles correspond to larger fiux densities. Note that a different size range for each channel was 
necessary for visualization purposes. 


small number of remaining sources, we do not create a PCCS2E 
for the LEI bands; instead, we flag the least reliable sources. 
Eurther details can be found in Sect. 3.2.1, which describes our 
assessment of the reliability of the PCCS2 at the lower frequency 
channels. 

HFI: There are no external full-sky catalogues at HEI frequen¬ 
cies; hence the reliability cannot be evaluated following the same 


procedure as for LEI. Instead, we follow a similar procedure to 
the one used for the PCCS as described in Planck Collaboration 
XXVIII (2014). In that paper, we used two measures of reliabil¬ 
ity for the HEI catalogues. The first measure, which we called 
simulation reliability, is determined from source injection into 
simulated maps and is defined as the fraction of detected sources 
that matched the positions of injected sources. If the simulations 
are accurate, such that the spurious and real detection number 
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Fig. 3. Distribution of the polarized sources in the lowest channels of the PCCS2. Red, green, and blue circles show sources from 
the 30, 44, and 70 GHz catalogues, respectively. As in Fig. 2, the size of the circle gives a qualitative idea of the relative polarized 
flux density of the source. 


counts mirror the real catalogue, this reliability is exact. The sec¬ 
ond measure, which we called injection reliability, is determined 
from source injection into the real maps, and in this case we 
are only simulating the real source component. The total source 
counts are composed of real and spurious components. In order 
to understand the reliability, we need to understand the spurious 
component. However, if we have a knowledge of the real compo¬ 
nent we can infer the spurious component from the total source 
counts. This is the motivation for this second approach, where 
injection reliability is deflned as the fraction of total extracted 
sources that correspond to injected sources recovered by PCCS2 
detection methods. This method therefore provides an approxi¬ 
mate assessment of the reliability. Naturally, the first method is 
preferred over the second, but the second approach becomes nec¬ 
essary when the simulations are not a sufficiently accurate repre¬ 
sentation of the sky signal. In order to establish whether the sim¬ 
ulations are of sufficient quality, we require that the simulated 
catalogues pass two internal consistency tests: that they have 
the same injected source completeness as the real catalogues; 
and that they have total detected number counts, as a function 
of S/N, that are consistent with those in the real data. (The in¬ 
trinsic number counts are assumed to be power law functions 
of flux density and are fitted to the detection counts at higher 
flux densities, where the catalogues are reliable and complete, 
and extrapolated to lower flux densities). In order to achieve this 
internal consistency we need to exclude the Galactic plane re¬ 
gion from the analysis, due to discrepancies arising from defi¬ 
ciencies in the simulation of diffuse dust emission near the beam 
scale (e.g. Galactic cirrus) and uncertainties in defining an input 
source model for the Galactic sources. The region excluded in¬ 
creases with frequency and is not a simple Galactic latitude cut, 
but is based on the level of the dust emission. The Galactic masks 
used at each of the six HFI frequencies are described in Sect. 3.2. 


In addition to these Galactic plane regions, for the highest four 
frequency channels we also exclude the region of sky inside a 
“filament mask” from the reliability assessment. Note that there 
is a different filament mask for each of these channels. These 
filament masks describe the areas of the sky in which residual 
structures, not related to sources, are present in the MHW2 Al¬ 
tered maps. The creation of these masks is explained further in 
Sect. 3.2.2. The union of the filament mask and the Galactic 
plane region then defines the area of the sky in which sources 
are assigned to the PCCS2E. Thus, whether a source is assigned 
to the PCCS2 or PCCS2E is determined solely by its location 
on the sky. The HFI reliability assessment for the PCCS2 is de¬ 
scribed further in Sect. 3.2.3 (recall there is no reliability assess¬ 
ment for the PCCS2E). 

2.4. Photometry 

As in the PCCS, we provide four different measures of the 
flux density for each source. They are determined by the 
source detection algorithm (DETFLUX), aperture photome¬ 
try (APERFLUX), point spread function fitting (PSFFLUX), 
and Gaussian fitting (GAUFLUX). Only the first is obtained 
from the Altered maps; the other measures are estimated 
from the full-sky maps at the positions of the sources. The 
source detection algorithm photometry, the aperture photome¬ 
try, and the point spread function (PSF) fitting use the Planck 
band-average effective beams, calculated with FEBeCoP (Fast 
Effective Beam Convolution in Pixel space; Mitra et al. 2011; 
Planck Collaboration IV 2016; Planck Collaboration VII 2016). 
Note that only the PSF fitting algorithm takes into account the 
variation of the PSF with position on the sky. The PCCS2 has 
been produced from the Planck full-mission maps (eight sky sur¬ 
veys in the LEI and five sky surveys in the HFI), and therefore su- 
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Fig. 4. Distribution of the polarized sources from the PCCS2 (top) and the PCCS2E (bottom). Red, blue, green, and black circles 
show sources from the 100, 143, 217, and 353 GHz catalogues, respectively. As in previous figures, the size of the filled circles 
gives a qualitative idea of the relative fiux densities of the sources at each frequency, where the larger dots correspond to larger fiux 
densities. Note that a different circle size range for each channel was necessary for visualization purposes. 


persedes the previous catalogues (for the PCCS only 1.5 surveys 
were analysed). It also includes the latest calibration and beam 
information, and we have improved some of the algorithms used 
to measure the photometry of the sources. In order to assess the 
differences between the photometry in the PCCS and PCCS2 we 
have compared both sets of catalogues at all Planck bands. 

A major change is in the Gaussian fitting photometry. We 
have implemented a new version of the algorithm that produces 
more robust measures, particularly for extended objects where 


the difference between the fiux densities in the PCCS and PCCS2 
can be as large as 100%. In the previous version of the al¬ 
gorithm, for some sources the fitting code was not converging 
properly and this issue has been addressed by using a new fit¬ 
ting approach; see Appendix B for a description of the method 
and its validation. In addition, the photometry from the detec¬ 
tion pipeline has changed at some frequencies by several per¬ 
cent, because it now takes into account the latest information 
about the effective beam FWHM and corrects for the biases 
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listed in Table 6, which range from 1% to 12%, depending on 
the frequency. The other two techniques, aperture photometry 
and PSF fitting, produce similar results in both catalogues. In 
the first case, the algorithm has not changed, while in the sec¬ 
ond, although the algorithm has been changed to improve the 
positional accuracy, this does not affect our measurements be¬ 
cause we use the coordinates from the detection pipeline as the 
reference for all photometric measures. In both cases the dif¬ 
ferences are always at the percent level. Moreover, flux densi¬ 
ties extracted from the publicly released Planck maps at 30, 44, 
and 70 GHz require a small correction for beam efficiency, since 
a small amount of power lies outside the main beam (Planck 
Collaboration II 2016). These small multiplicative corrections 
are 1.00808, 1.00117, and 1.00646 at 30, 44, and 70 GHz, re¬ 
spectively. The flux densities provided for sources in the PCCS2 
catalogues have been corrected accordingly. Uncertainties are 
provided for all four flux-density measures. In Table 13 we show 
the uncertainties associated with the flux densities of the faintest 
sources in the extragalactic zone of each catalogue, after exclud¬ 
ing the faintest 10 % of sources as obtained with the source de¬ 
tection algorithm. These uncertainties range from 90 to 130mJy 
for the 30-70 GHz catalogues, and from 30 to 270 mJy for the 
100-857 GHz catalogues. This gives an idea of the sensitivity 
of the catalogue and the associated uncertainties. However, the 
uncertainties depend not only on the fiux density of the sources 
but also on their position in the sky, so we provide noise maps 
that can be used to estimate the expected uncertainty in the fiux 
density of a source at any position in the sky. 


Detection pipeline photometry (DETFLUX). The detection 
pipelines assume that sources are point-like. The amplitude of a 
detected source is converted to fiux density using the solid angle 
of the effective beam (from Table 2), and the conversion from 
map units into intensity units. The uncertainty in the fiux den¬ 
sity for each source is measured as the local noise in an annulus 
around the source in the MHW2-filtered map, where bright pix¬ 
els belonging to other compact sources in the vicinity, if any, are 
excluded from the calculation. If a source is resolved by Planck 
its fiux density will be underestimated. In this case it may be 
better to use the GAUFLUX estimation. The estimation of the 
fiux density provided by the HFI detection pipeline has been im¬ 
proved since the PCCS release, by removing a bias that lowered 
the recovered fiux densities in the higher frequency channels (see 
Planck Collaboration XXVIII 2014). The photometric perfor¬ 
mance of the PCCS2 detection pipeline is assessed in Sect. 3.4. 


Aperture photometry (APERFLUX). The fiux density is esti¬ 
mated by integrating the data in a circular aperture centred at 
the position of the source. An annulus around the aperture is 
used to evaluate the level of the background. The annulus is also 
used to make a local estimate of the noise to calculate the un¬ 
certainty in the fiux density. The fiux density is corrected for the 
fraction of the beam solid angle falling outside the aperture and 
for the fraction of the beam solid angle falling in the annulus. 
The aperture photometry was computed using an aperture with 
a radius equal to the average FWHM of the effective beam (the 
effective FWHM in Table 2), and an annulus with an inner ra¬ 
dius of 1 FWHM and an outer radius of 2 FWHM. The effective 
beams, also given in Table 2, were used to compute the beam 
solid angle corrections. For details see the PCCS paper (Planck 
Collaboration XXVIII 2014). 


PSF fit photometry (PSFFLUX). The fiux density and its un¬ 
certainty are obtained by fitting a model of the PSF at the po¬ 
sition of the source. The model has four free parameters: the 
amplitude of the source; a background offset; and two coordi¬ 
nates for the location of the source. The PSF is obtained from 
the effective beam by means of a bicubic spline interpolation for 
source positions that are different from the centre of a pixel. Note 
that the PSF fitting now includes subpixel positioning, which is 
a new feature introduced after the production of the PCCS. For 
details see Appendix A. 

Gaussian fit photometry (GAUFLUX). The approach to 
Gaussian fitting has been completely revised since the PCCS. 
The algorithm now allows the position of the source to vary 
as the best fit is found. The same parameters are returned for 
each source: its flux density; the major and minor semi-axes; 
and an orientation angle. Additionally, as in the PCCS, the semi¬ 
axis values are used in the construction of the flag for extended 
sources. The new method uses a downhill simplex method in 
multidimensions, the Nelder-Mead method, to find the best-fit 
values in the full parameter space of position, fiux density, and 
elliptical Gaussian parameters. The method has been shown to 
be robust and stable (Press et al. 1992). Optimization is based on 
the reduced log-likelihood with prior regularization for the size 
of the source defined by the effective beam at each frequency. 
The downhill simplex methods does not produce estimates of 
the fiux density uncertainties. For this purpose a Markov Chain 
Monte Carlo has been used. For details see Appendix B. 

2.5. Polarization 

In the Planck polarization maps, the polarized sources are em¬ 
bedded in a background that is the combination of instrumen¬ 
tal noise and diffuse emission. The nature of the diffuse emis¬ 
sion depends on the observation frequency; for example, po¬ 
larized synchrotron emission in the lower frequency channels 
and infrared emission in the higher frequency channels. In both 
regimes the polarization fraction of the compact sources (the ra¬ 
tio between their polarized fiux densities and total intensity) is 
typically lower than 1-2%. This presents a challenge in terms 
of disentangling the true polarized fiux density of a source from 
the background. In order to tackle this problem, a two-step pro¬ 
cess has been proposed (Lopez-Caniego et al. 2009). First, a 
maximum-likelihood filter is applied, reducing the noise and en¬ 
hancing the S/N of the sources embedded in the Q and U maps 
(Argueso et al. 2009). Second, the significance of each detec¬ 
tion is assessed based on the statistics of the local background 
in the vicinity of the source. Several significance levels were in¬ 
vestigated and we concluded that, for the typical polarization 
backgrounds present in the Planck polarization maps, a signif¬ 
icance threshold of 99.99 % successfully distinguishes the po¬ 
larized emission of a compact source from a peak in the back¬ 
ground. This approach has been used in the present catalogues to 
attempt to measure the polarized fiux densities and uncertainties 
of all sources found in the temperature maps. Polarization mea¬ 
surements are provided for all sources where the significance 
of the detected polarized signal reaches or exceeds the limit of 
99.99 %; for the remaining sources we provide the 99 % upper 
limit. Figures 3 and 4 show the distributions of the significantly 
polarized sources in the LFI and HFI polarized frequency chan¬ 
nels. Normalized histograms of the polarization fraction for the 
population of significantly polarized sources in the PCCS2 cata¬ 
logues are shown in Fig. 5. 


7 


Planck Collaboration: Second Planck Catalogue of Compact Sources 



0 5 10 15 20 


Polarization fraction [%] 


Fig. 5. Normalized histograms of the recovered polarization 
fraction from the PCCS2 catalogue from 30 to 217 GHz (the 
353 GHz channel is not shown because the catalogue contains 
only one source). The number of sources in each histogram is 
indicated. 


If a source is not strictly point-like, the filtering proce¬ 
dure used to reduce the noise will also remove signal. For this 
reason we also provide aperture photometry measurements for 
polarization, which, although noisier, do not remove as much 
signal from compact (but not point-like) sources as filtering. 
The aperture photometry package is common to both LFI and 
HFI, whereas there are two different implementations of the 
maximum-likelihood estimator, which will allow us to assess the 
robustness of the methods by comparing their results in a com¬ 


mon set of simulations. The results of the HFI and LFI polariza¬ 
tion pipelines are compared in Sect. 3.5. 

The polarized flux density of a source, F, is evaluated using 

p = (1) 

where Q and U are the flux densities in the Stokes Q and U 
maps, measured at the position of the source detected in the 1 
map. We follow the lAU/IEEE convention (Hamaker & Bregman 
1996) for defining the angle of polarization of a source: polariza¬ 
tion angles are taken as increasing anticlockwise (north through 
east). In this paper, position angle zero is taken as the direction 
of the north Galactic pole. The polarization angle is defined by 

i//= iarctan(-f//e). (2) 

The minus sign is necessary to correct from the HEALPix con¬ 
vention for position angles used in the Planck Stokes parameter 
maps, in which position angle increases clockwise. As noted in 
Planck Collaboration I (2016), the convention used for Planck 
polarization maps is the one usual in CMB studies and is used 
in WMAP papers, whereas the lAU/IEEE convention adopted 
here is standard for astronomical sources. Polarization angles are 
given in degrees in the range -90° to 90°. The estimate of P ac¬ 
quired using Eq. (1) is biased, because the errors in the Q and 
U measurements, on average, contribute positively to the mea¬ 
surement of P (see, e.g. Montier et al. 2015). However, in our 
significance regime we can use the approximation 

^debiased — -yj~ (5) 

to debias our estimate of P, where is the error in P and is 
calculated by propagating the errors in Q and U, where ctq^u 
are calculated as the local rms in an annulus around the source 
in the maximum-likelihood filtered Q and U maps, under the 
assumption of no correlation: 



The polarization angle error is obtained by propagating the er¬ 
rors in Q and U : 

" 2(6^ +1/2) ^ 

As shown in Table 14 the typical uncertainty in the polar¬ 
ized flux density is 45-90 mJy between 30 and 70 GHz and 30- 
ISOmJy between 100 and 353 GHz. 

2.5.1. Corrections for bandpass mismatch 

Mismatch between the bandpass shapes of the two orthogonally- 
polarized detectors in each feed horn causes leakage of total in¬ 
tensity into the polarization signal for any emission whose spec¬ 
trum differs from that of the primary calibrator, namely the CMB 
dipole; therefore all foreground emission including that from 
compact sources suffers from temperature-to-polarization leak¬ 
age. Correction requires a model of the spectrum of the source, 
as well as a model for the spectral response of each detector 
or bolometer. Since the detecting elements used in the two in¬ 
struments are different, LFI and HFI treated bandpass mismatch 
differently. The magnitude of the correction can be very differ¬ 
ent from one source to another. In the lower Planck frequencies, 
the correction can vary from a fraction of a percent up to 100 %. 
In the higher frequency channels this correction is always below 
the percent level. The details are presented in Appendix C. 
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2.5.2. Evaluation of marginal polarization measurements 

At four HFI frequencies (100, 143, 217, and 353 GHz) we 
present an additional set of polarized flux-density and polariza¬ 
tion angle estimates for sources detected only marginally in po¬ 
larization. These are derived using the Bayesian “PowellSnakes” 
algorithm (see Planck Collaboration XXVIII 2014). The aim is 
to disentangle the sources that have some polarized emission 
from those that are consistent with no polarized signal. This al¬ 
lows us to probe fainter polarization signals, and thus to provide 
deeper and more complete polarization catalogues without any 
loss of reliability (as we show in Sect. 3.5). The details of the 
method are presented in Appendix C.3. 

3. Validation of the PCCS2 

The contents of the PCCS2 and the four different flux-density es¬ 
timates have been validated by simulations (internal validation) 
and comparison with other astrophysical data (external valida¬ 
tion), as was done for the PCCS (Planck Collaboration XXVIII 
2014). The validation of the low-frequency sources can be per¬ 
formed in part by using the large number of existing catalogues. 
Detections identified with known sources have been flagged as 
such in the catalogues. In contrast, the validation of sources at 
higher frequencies must be done using simulations, specifically 
through a Monte Carlo quality assessment process in which arti¬ 
ficial sources are injected into both real and simulated maps. In 
the following subsections, we discuss tests on the completeness, 
reliability, astrometry, and photometry of the single-frequency 
catalogues, as well as comparisons between different Planck 
bands. We also describe internal and external validation of the 
polarization measurements. 

3.1. Completeness 

3.1.1. LFI 

In the case of the three lowest frequencies, we compared the 
Planck compact source detections to external catalogues of ra¬ 
dio sources. We began by constructing a band-merged catalogue 
based on positional coincidence and included all the sources 
detected above the initial 4(T detection threshold at any of the 
LFI frequencies. The catalogue contained 2039 sources, many 
of which were detected at only one frequency. We then used 
this catalogue to make a position-based search for identifica¬ 
tions with three external catalogues of radio sources using a 
search radius of 1.5 x (Xb, where (Xb = FWHM/(2 V2 In 2) can 
be evaluated from the fitted FWHM in Table 2. These catalogues 
are: (1) in the southern hemisphere, the Australia Telescope 
20 GHz Survey (AT20G; Murphy et al. 2010), a catalogue of 
sources brighter than 40 mJy that covers the whole southern sky 
(\b\ < 0°); (2) in the northern hemisphere, where no large-area 
survey at similar frequencies to AT20G is available, the 8.4 GHz 
Combined Radio All-sky Targeted Eight GHz Survey (CRATES; 
Healey et al. 2007), a compilation of flat-spectrum (a> - 0.5) 
radio sources with nearly uniform extragalactic (\h\ > 10°) 
coverage for sources brighter than 65 mJy at 4.8 GHz; (3) the 
full-sky New WMAP Point Source Catalogue (NEWPS; Lopez- 
Caniego et al. 2007; Massardi et al. 2009) covering the frequency 
range 23-61 GHz, that include sources brighter than 700 mJy at 
23 GHz and is complete above 2 Jy. These catalogues have a sim¬ 
ilar source density, so there should not be multiple associations 
of AT20G, CRATES, or NEWPS sources within a Planck beam, 
and their frequencies range from 8 to 61 GHz. The NEWPS cat- 



Flux density [Jy] 



Fig. 6. Validation of the sources in the 30, 44, and 70 GHz chan¬ 
nels. The top panel shows the cumulative completeness per flux- 
density bin of the PCCS2 catalogue. The lower panel shows 
the number of unidentified sources per flux-density bin. These 
sources have no clear counterpart in any of the external cata¬ 
logues used for validation. However, many of these source have 
been detected in more than one Planck channel and this implies, 
first, that they could be real sources and not spurious, and one 
can consider the number of non-matched sources presented here 
as an upper limit. Second, the fact that they are not in the ex¬ 
ternal catalogues suggests that, if real, they could be potentially 
interesting objects, maybe going through a flaring phase. 


alogue was produced by analysing the WMAP maps, and this 
dataset is very similar to that of Planck in terms of format (all¬ 
sky Healpix maps), background characteristics, and angular res¬ 
olution (-13-33' in Planck LFI, vs. -14-56' in WMAP), which 
makes it a very good dataset for validation. In addition, the two 
catalogues, PCCS2 and NEWPS, were produced using the same 
source-extraction tool (the Mexican Hat wavelet). For the search 
radius, we use 1.5 x (Xb, which includes -87 % of the area of the 
Planck beam, as a compromise between the arcsecond resolu- 
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tion of CRATES and AT20G, and the ~ degree resolution of the 
23 GHz WMAP channel. In any case, given the similar source 
density of the catalogues, one could use a slightly larger search 
radius and the results would not change. As in the PCCS, the 
PCCS2 includes 94% of the sources in NEWPS when using a 
1.5 X (T^ search radius, so in order to study the completeness 
of the catalogue deeper samples like CRATES and AT20G are 
needed. However, the frequencies of these two surveys, 8.4 and 
20 GHz, are lower than the lowest Planck frequency, and vari¬ 
ability and spectral effects could push some of the sources below 
the PCCS2 detection thresholds. Thus, the completeness that we 
estimate by comparing the PCCS2 against these three catalogues 
is a lower limit. Eor this reason we used an alternative complete¬ 
ness estimate that can be derived from knowledge of the noise 
in the maps when the completeness is greater than 50%. If the 
flux density estimates S are subject to Gaussian errors with am¬ 
plitude given by the noise of the Altered patches, the cumulative 
completeness per patch should be 


C(5)=I + Ierf 


S - go-si0,(k) \ 
crs(0,(p) /’ 


( 6 ) 


where crj ( 6 , (p) is the variance of the Altered patch located at 
( 6 , (p),q is the S/N threshold and erf(x) = -^ dt is the stan¬ 
dard error function. The true completeness will depart from this 
limit when the simplifying assumptions of non-Gaussian noise 
and uniform Gaussian beams are broken. Using this expression, 
the cumulative completeness of each LEI band is derived by 
making use of a model of the source counts N{S) (de Zotti 
et al. 2005) that accounts for various source populations (flat- 
spectrum radio quasars, BE Lac objects, steep-spectrum sources, 
GPS sources, early phase gamma-ray after glows, etc.). The 
true completeness will depart from this limit when the simpli¬ 
fying assumptions of non-Gaussian noise and uniform Gaussian 
beams are broken. 

Eigure 6 demonstrates that the catalogues at 30, 44, and 
70 GHz are essentially 100% complete above flux densities of 
1 Jy, and for 30 and 70 GHz are still more than 95 % complete 
down to around 0.6 Jy. The flux densities that correspond to the 
90 % completeness level are shown in Table 13. 


3.1.2. HFI 

There are no external full-sky catalogues in the HEI frequency 
range, so we rely on the injection of artiflcial sources into the 
Planck maps to establish the completeness, following the power- 
law models fitted to the data described in Sect. 3.2.3 and Table 4. 

The completeness is determined from the injection of unre¬ 
solved point sources into the real maps. Bias due to the super¬ 
imposition of sources is avoided by preventing injection within 
an exclusion radius of ct/,, evaluated from the fitted EWHM in 
Table 2, around both existing detections in the real map and pre¬ 
viously injected sources. We note that while superimposition of 
real sources will occur in the Planck maps, it is beyond the scope 
of this paper to disentangle such effects. Our definition of a point 
source here is a beam-shaped spike of emission, regardless of the 
make-up of astrophysical objects that produce the emission. 

The flux from real and injected point sources contributes to 
the noise estimation for each patch, reducing the S/N of all de¬ 
tections and biasing the completeness. We prevent this effect by 
determining the noise properties on the maps before injecting 
sources, and have verified that any remaining bias on detection 
and parameter estimates due to injected sources is negligible. 


The injected sources are convolved with the effective beam com¬ 
puted using the FEBeCoP algorithm (Mitra et al. 2011; Planck 
Collaboration IV 2016; Planck Collaboration VII 2016). 

We show the Monte Carlo completeness functions computed 
per channel in Pig. 7. We have calculated them using several 
cuts based on the “highest reliability catalogue” column that in¬ 
dicates the highest reliability catalogue to which each source be¬ 
longs (see Sect. 3.2.3); as the required catalogue reliability in¬ 
creases, the completeness decreases. We also show the complete¬ 
ness for the unvalidated area of the sky containing the PCCS2E. 
This area is significantly less complete than the main PCCS2 
regions due to Galactic emission, where the completeness dete¬ 
riorates markedly with frequency. The PCCS2 product includes 
noise maps at each channel, which are also shown in the mid¬ 
dle panels of Pig. 8. These maps contain the MHW2 detection 
noise, defined as the standard deviation of the MHW2 Altered 
patches, 0 -^( 6 , (p). A good approximation to the Monte Carlo 
completeness can be calculated for any subsection of the sky us¬ 
ing the noise maps and the reliability threshold maps. Assuming 
Gaussian noise in the Altered patches, the completeness is given 
by Eq. (6). We compare the error function and Monte Carlo com¬ 
pleteness estimates for the 80 % reliability catalogues at each 
channel in Pig. 9. The effective that we used, and included 
in the data release, has been normalized across the PCCS2 re¬ 
gion to match the effective noise from the Monte Carlo tests. 
The Monte Carlo completeness drop-off is slightly wider than 
the error-function completeness because it includes the effects 
of non-Gaussian noise from the background and varying asym¬ 
metric beams. However, the discrepancy is less than 5 %. 

3.2. Reliability 

The underlying philosophy of the PCCS2 was to provide a cata¬ 
logue of sources with a reliability of at least 80 %, which could 
be cut if desired by a user to generate higher reliability subsets. 
This has required, as explained in Sect. 2.2, the division of the 
HEI data into two subcatalogues, the PCCS2 and PCCS2E. The 
assignment to the PCCS2E is based on the position of the source 
on the sky. These locations are determined by the union of the 
Galactic region and the filament mask. Note that at 857 GHz this 
union represents more than half of the sky. The PCCS2E there¬ 
fore contains substantially more sources than the PCCS2 in the 
higher frequency channels since it includes all the Galactic plane 
sources. 


3.2.1. LFI: reliability assessment 

The band-merged catalogue compared with ancillary data, as de¬ 
scribed in Sect. 3.1.1, has also been used to assess the reliabil¬ 
ity of the LFI catalogues. All but 426 PCCS2 sources (out of 
2039) were identified with known sources in the external ra¬ 
dio source catalogues of CRATES, NEWPS, and AT20G. The 
initial estimate for the reliability, using just these external cata¬ 
logues, is thus greater than 79 %. The percentage of unidentified 
LFI sources as a function of their flux density is shown in the 
lower panel of Fig. 6. Associations for some of these initially 
unidentified sources were later found by performing a search 
in NED. This procedure was carried out on a source-by-source 
basis and some subjective judgement was required. Hence, we 
report a range in the possible number of sources thus identi¬ 
fied. The sources that are positively identified by this approach 
are flagged as such and the associated sources are named in the 
PCCS2 catalogue. 
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Fig. 7. HFI completeness results from the Monte Carlo quality assessment. Completeness is shown per channel for the 80, 85, 90, 
and 95 % reliability catalogues, and for the PCCS2E. 


Reliability at high Galactic latitudes using external cata¬ 
logues: out of the 426 initially unidentified sources, 180 were 
at high Galactic latitudes, \h\ > 20°. More than a quarter of these 
were positively identified using NED, leaving 132. Given the 
1161 sources detected over this latitude range, this implies a re¬ 
liability of 92 % for high Galactic latitude sources. It should be 
noted that 7-8 % of the remaining 132 unidentified sources ap¬ 


pear in multiple Planck bands, and are therefore likely to be real. 
Hence, this reliability estimate is a lower bound. 

Reliability using external catalogues and taking into account 
multiple band detections: we also searched for identifications for 
all sources that appeared in two or more Planck bands; here we 
included the 100 GHz band in the analysis. Of the 426 initially 
unidentified sources, 133 appeared at two or more frequencies. 
Of these, the number that remained unidentified was 71. Given 
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Fig. 8. Left: zone masks constructed from the filament masks and the Galactic masks. The area covered by the PCCS2 is given by 
the zero (blue) values in the mask, and that covered by the PCCS2E is given by the non-zero values. The Galactic region is traced 
in red, and the filament mask is yellow outside the Galactic region and green inside it. Centre: rms noise level as determined by the 
HFI MHW2 code. Right: S/N thresholds applied to the raw catalogue; a fiat S/N cut is applied in the region of the PCCS2E and for 
the PCCS2 at 545 and 857 GHz. 
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Input flux density (Jy) 

Fig. 9. Comparison of the error-function semi-analytic com¬ 
pleteness derived from the threshold maps (Eq. 6) and the 
Monte Carlo completeness estimates. Top: solid lines denote 
the error-function completeness and the dashed lines the Monte 
Carlo completeness. Bottom: difference between the two com¬ 
pleteness estimates. 


that more than half of the LFI sources (1149) appear at two or 
more frequencies, this implies a reliability > 94 % for sources 
detected in two or more Planck bands. 

Reliability using external catalogues and sources detected in 
a single band and with SNR > 5: we searched NED for identi¬ 
fications for sources detected in only a single Planck band with 
S/N > 5. Not surprisingly, the rate of identification was lower, 
but a few positive identifications were added. In the end we were 
left with 335 out of the original 426 that remained either un¬ 
examined or unidentified. We may therefore conclude that the 
overall reliability of the LEI catalogue is at least 84 %. 

The procedure described above was used to construct the 
PCCS2 for the three lowest Planck bands, and additionally to 
populate the EXT_VAL column for them. This column summa¬ 
rizes the comparison with external catalogues, and is described 
in Sect. 5. Note that the higher-frequency bands also provide in¬ 
formation for each source in this column, but they are not used 
in the reliability assessment. 


3.2.2. HFI: Filament masks 

There are regions on the sky even at high Galactic latitudes in 
which the detection of sources cannot be trusted to high levels 
of reliability. The reason is that the Mexican hat wavelet algo¬ 
rithm used for the detection of compact sources is also efficient 
at edge detection. While this is ideal for point sources embed¬ 
ded in a Gaussian noise background, it is not optimal for a map 
with non-Gaussian structures where there are other edges to de¬ 
tect, such as dusty filaments. The impact of false detections due 
to these structures may in part be ameliorated by rejection crite¬ 
ria based on the number of connected pixels of the detection, as 
described in Planck Collaboration XXVIII (2014); however, this 
is not a complete solution. We wish to place source detections 
occurring in these filamentary structures into the PCCS2E rather 
than the main PCCS2. In order to do this we need to create fil¬ 
ament masks, which describe the regions of the sky containing 
filamentary structures that pass into the wavelet filtered patches. 

The obvious way to construct these masks would be to use 
the filtered patches themselves, after median filtering to remove 
the point-like objects from the patch. The difficulty with this ap¬ 
proach is setting an appropriate threshold above which to mask. 
Given that these dusty structures are far from uniformly dis¬ 
tributed across the sky, local evaluations of this threshold will 
not be as successful as a global one. Any practical filament mask 
should also be continuous; for both these reasons they cannot be 
created directly from the filtered patches. 

A Mexican hat wavelet filtered map is very close to a 
difference-of-Gaussians (DOG) map if the ratio of the two 
EWHMs used is 2:1. The DOG map is created by smoothing 
the original map with two different Gaussians, creating two new 
maps that are then differenced.^ Here the smaller of the two 
EWHMs used is the fitted FWHM from Table 2. A full-sky DOG 
map may be trivially created, and while it is not identical to 
the filtered patches it traces the same structures in the original 
map that pass into the filtered patches. In order to remove point 
sources, the DOG map is median-filtered using a filter radius 
of twice the fitted FWHM (from Table 2). By thresholding this 
median-filtered DOG (MF-DOG) map, we can create the desired 
filament mask for each channel. In order to find an appropriate 
threshold, we select the cleanest (faintest) 25 % of the sky based 
on the smoothed sky brightness of the 857 GHz channel map, 
and create a histogram of the MF-DOG map from the pixels in 
this region. We then fit a Gaussian to this histogram, and the 
filament mask is given by all the pixels in the MF-DOG map 
with values greater than 3 times the a of this fitted Gaussian. 
Negative fluctuations are not masked, since they cannot lead to 
spurious detections. In the case of point sources in a Gaussian 
background, this procedure would result in an MF-DOG map 
that would contain solely Gaussian noise, and a mask created 
as above would be expected to mask 0.15 % of the sky. Table 3 
shows the percentage area of the filament mask for the frequency 
channels at which it is used; we see that the percentage area 
masked is well in excess of the expectation from Gaussian statis¬ 
tics. Additionally, we see that the area of the filament mask in¬ 
creases with frequency as the maps contain more and more emis¬ 
sion from dust. Indeed, we do not need to use a filament mask 
for the lower two HFI frequency channels, because the dusty fil¬ 
amentary structures are not a problem at these frequencies. 

Note that the risk of high S/N point sources located inside 
filamentary structures being placed inside the filament mask is 
limited by the fact that high S/N point sources in the DOG map 

^ This is the default approach used for finding compact sources in 
data from the SCUBA-2 instrument (e.g. MacKenzie et al. 201 1). 
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Table 3. Fractions of sky covered by the PCCS2 catalogue, the 
filament mask, and the Galactic region for the HFI channels. No 
masks are required for the LFI channels, where the PCCS2 cov¬ 
ers the entire sky. The filament mask is not required for the low¬ 
est two HFI frequency channels. Some parts of the sky are cov¬ 
ered by both the filament mask and the Galactic region. 


Channel 

PCCS2 area 
% sky 

Filament mask 
% sky 

Galactic region 
% sky 

100 ... . 

85.0 


15.0 

143 ... . 

85.0 


15.0 

217 ... . 

64.9 

2.2 

35.0 

353 ... . 

47.6 

7.5 

52.0 

545 ... . 

47.0 

10.9 

52.0 

857 ... . 

46.3 

13.9 

52.0 


Table 4. Power law model fit parameters obtained from the sim¬ 
ulation of the source number counts in the HFI channels used in 
the reliability assessment. 


Channel 

*^min 

[Jy] 

taper 

[Jy] 

a 

100 ... . 

0.3 

0.1 

2.54 

143 ... . 

0.3 

0.2 

2.51 

211 ... . 

0.2 

0.1 

2.63 

353 ... . 

0.4 

0.1 

2.69 

545 ... . 

0.7 

0.1 

2.59 

857 ... . 

1.5 

0.3 

2.34 


are positive peaks surrounded by negative troughs. This means 
that once the DOG map is median-filtered the resultant level in 
the MF-DOG map at the location of these point sources is likely 
to be below the level at which thresholding occurs. 

3.2.3. HFI: reliability assessment 

As described in Sect. 2.3, we use two methods to assess the HFI 
reliability, the simulation reliability and the injection reliability. 
The simulation reliability is assessed by injecting sources into 
simulated maps; this simulates both the real and spurious de¬ 
tection components of the total counts. The injection reliability, 
however, involves injecting sources into the real maps; hence 
only the real detection component of the total counts is simu¬ 
lated. The real component to inject, in both cases, is evaluated 
using the PCCS2 above a given fiux density threshold (above 
which we are complete) fitted to a single power-law model de¬ 
fined as N{S) oc S~^. For each channel we provide three num¬ 
bers in Table 4: a is the power-law index (estimated per fre¬ 
quency with error in the range 0.05-0.1); 5 min is the minimum 
fiux density considered when fitting the model (i.e. lower fiux 
densities are excluded to allow for incompleteness); and 5taper 
is the fiux density at which the power law was truncated for 
the completeness simulations to avoid dominating the injected 
population with unobservable faint sources. Note that there was 
no truncation for the reliability work (simulation or injection). 
Simulation reliability is preferable to injection reliability be¬ 
cause it provides a more complete understanding of the detec¬ 
tion properties of the catalogue, including information on how 
the reliability varies as a function position on the sky, as well as 
on S/N. 

Previously, for the PCCS (as described in Planck 
Collaboration XXVIII 2014), the simulation reliability was used 


for 100-217 GHz, while injection reliability was used for 353- 
857 GHz. This was because the simulated maps at 353-857 GHz 
could not be used to produce a simulation reliability estimate. 
The simulated number counts and completeness were not con¬ 
sistent with the real data, and the discrepancies arise from de¬ 
ficiencies in the simulation of diffuse dust emission near the 
beam scale. Since then, improvements have been made in the 
simulations; using the FFP8 simulations (Planck Collaboration 
XII 2016), it is now possible to extend the use of simulation 
reliability to 353 GHz. These maps, however, include a leaked 
compact-source component from the Planck maps from which 
they were derived. As this can produce artificial, high S/N, spu¬ 
rious sources, we screen these from the reliability estimates by 
considering any detection at S/N > 10 to be real. At 545 and 
857 GHz, we continue to use the injection reliability estimate. 

Table 3 shows the percentage area of the sky occupied by the 
PCCS2, for each HFI channel. This corresponds to the area in 
which the reliability assessment is performed. The PCCS2 cov¬ 
ers the region of the sky not excised by the filament mask or 
the Galactic region. The Galactic region is determined by the 
area, from each channel, that must be excluded from the relia¬ 
bility assessment in order to achieve consistency between simu¬ 
lated and real catalogues. The percentage areas of the Galactic 
regions and filament masks are also shown in Table 3. Note that 
the Galactic region and the filament mask can and do overlap; 
their union forms the area of sky in which sources are assigned 
to the PCCS2E. Hence, the complement of the PCCS2 is the 
PCCS2E, and consequently no source can appear in both sub¬ 
catalogues. 

Each entry in the PCCS2 and PCCS2E catalogues contains a 
field, HIGHEST_RELIABILITY_CAT, which specifies the high¬ 
est reliability subsample in which that source may be included. 
Hence, this field can be used to cut the catalogues to produce 
subsets with a reliability higher than the survey target of 80%. 

Eor the 100-353 GHz channels we can perform the reliabil¬ 
ity assessment using simulation reliability. This allows us to de¬ 
fine a local S/N threshold, q{6, 0, !P), as a function of both sky 
position (jd, 0) and target reliability, P. This threshold gives a lo¬ 
cal reliability P within an 8° radius of {6, 0). This information 
allows us to populate the HIGHEST_RELIABILITY_CAT field 
for all sources in the PCCS2 catalogues in these channels, and a 
resolution of 1 % reliability. Note that, for these four channels, 
the option to create a higher reliability subset will also apply to 
spatial subsets of the original catalogue. Eor example, one could 
create a catalogue of the north ecliptic pole region to a reliability 
of 97 %, if desired. 

Eor the 545 and 857 GHz channels the limitations of the 
injection reliability mean that we are unable to define a lo¬ 
cal S/N threshold for a given reliability. Hence, we provide 
a global S/N threshold that will deliver the target reliability 
for the full catalogue. We use this approach to populate the 
HIGHEST_RELIABILITY_CAT field for these channels in steps 
of 5 % in reliability. Since there is no local assessment, the op¬ 
tion to create higher-reliability spatial subsets for these channels 
is not available, because the desired reliability will only apply to 
the full area covered by the catalogue. We note that the PCCS2 
survey S/N threshold at 857 GHz is substantially higher than 
the threshold applied to build the PCCS. We have improved the 
modelling of the real extragalactic sources for the injection relia¬ 
bility estimate for the PCCS2, which results in a shallower spec¬ 
tral index for the input source model and fewer injected sources 
at low S/N. This produces a more realistic, and lower, reliability 
estimate at a given S/N relative to the PCCS. As the threshold 
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has moved to higher S/N, the flux density at 90 % completeness 
in Table 13 is now higher than for the PCCS. 

The S/N threshold maps which produce the survey target of 
80 % integral reliability are shown in the right-hand panels of 
Fig. 8, for all the HFI channels. Note the flat S/N cut applied 
to the 545 and 857 GHz channels in the PCCS2 regions, as well 
as the flat S/N = 5 cut applied to the PCCS2E regions. The S/N 
threshold HFI maps for !P = 80%, 85 %, 90%, and 95 % are 
included in the data release. 

3.2.4. HFI: Comparison with H-ATLAS 

As an external check on the reliability of the PCCS2 we have 
exploited the catalogue of submillimetre sources extracted from 
the full Herschel Astrophysical Terahertz Large Area Survey (H- 
ATLAS; Bales et al. 2010). The survey covers an area of about 
550 deg^ at 250, 350, and 500 /im (the two longer wavelength 
channels corresponding fairly closely with HFFs 857 GHz and 
545 GHz). A public catalogue is available only for the H-ATLAS 
Science Demonstration Phase held (covering ^ 16deg^) so we 
have used, with permission, one that is unpublished and was re¬ 
leased for internal use of the H-ATLAS consortium (Maddox et 
al. 2015; Valiante et al. 2015). Within the H-ATLAS flelds, the 
PCCS2 catalogue contains 39, 44, and 121 sources at 353, 545, 
857 GHz, respectively, while the PCCS2E contains two sources, 
both detected only at 857 GHz. We have identifled, at each fre¬ 
quency, the PCCS2 with the H-ATLAS catalogue using a search 
radius equal to half the Planck FWHM. Increasing the search ra¬ 
dius to one FWHM does not add any additional reliable counter¬ 
parts. As expected, given the large surface density of H-ATLAS 
sources, at least one source is always present within the search 
radius. We have taken as reliable counterparts to Planck sources 
those with H-ATLAS flux densities within a factor of three of the 
Planck APERFLUX ones. Since the Herschel photometry does 
not extend to 353 GHz, to look for reliable counterparts at this 
frequency we have extrapolated the H-ATLAS flux densities us¬ 
ing the spectral properties measured at higher frequencies. In the 
present, preliminary, version of the H-ATLAS catalogue, aper¬ 
ture corrections of flux densities for extended sources have been 
applied only for a fraction of the area. It is thus possible that a 
few H-ATLAS flux densities of large galaxies have been strongly 
underestimated and are therefore missed as reliable counterparts 
to Planck sources by the flux-density criterion. To recover them 
we have looked for associations of Planck sources with large, 
nearby optical galaxies. In this way we have recovered flve H- 
ATLAS counterparts at each of the Planck frequencies. We And 
26, 38, and 112 reliable counterparts to PCCS2 sources at 353, 
545, and 857 GHz, respectively, which translates into a reliabil¬ 
ity of 26/39 (67 %X 38/44 (86 %), and 112/121 (93 %). Neither 
of the two sources in the PCCS2E has a reliable H-ATLAS coun¬ 
terpart. 

The 353 GHz channel is flve identiflcations short of the target 
80 % integral reliability. Since we require consistent flux densi¬ 
ties between Planck and the extrapolated values from Herschel 
for an identiflcation, this slight deficit should not raise too much 
concern about the reliability of the 353 GHz catalogue. 

An assessment of the expected number of matches due to 
chance between H-ATLAS and Planck sources was made using 
the number of Planck detections within the H-ATLAS area, the 
size of the Planck beams, and the number density of H-ATLAS 
sources with flux densities above the 90% completeness lim¬ 
its of the PCCS2 (scaling these limits by a factor of three, be¬ 
cause we have taken as reliable counterparts to Planck sources 
those with H-ATLAS flux densities within a factor of three of 


the Planck APERFLUX ones). It was found that the number of 
random associations is < 10“"^ for all frequencies studied ( 353, 
545, and 857 GHz). 

3.3. Astrometry 

The astrometric accuracy and positional uncertainties of the 
PCCS2 and PCCS2E were determined using both internal and 
external tests. The external validation was based on a compar¬ 
ison of PCCS2 source positions with those measured with the 
Karl G. Jansky Very Large Array (VLA), which in turn are tied 
to the ICRF reference frame to an accuracy well below 1". This 
comparison both validates the astrometric accuracy of the Planck 
catalogues and provides estimates of positional uncertainties, 
channel by channel. We also used simulations, in which we in¬ 
jected sources to test the accuracy of the positions determined 
by the detection algorithms. 

3.3.1. External tests 

For the six lowest frequency channels, a direct comparison was 
made between the PCCS2 coordinates of bright, non-thermal ra¬ 
dio sources and the subarcsecond precision positions determined 
by the VLA. Table 5 summarizes the observed offsets in position 
between Planck and VLA observations for up to 70 (depend¬ 
ing on the frequency) compact, unconfused synchrotron sources. 
The positional offsets recorded in Table 5 are the averages (to¬ 
gether with the error on the mean) of all the sources in this study; 
hence these results limit any global positional offset between 
the VLA and Planck, and do not represent the positional off¬ 
set for a single source. The positional errors given in Table 5 
are evaluated from the standard deviations of the offsets in lat¬ 
itude and longitude, and hence provide the radial error in the 
position of an individual source in each channel. The numbers 
of sources used at each channel are also given. We also exam¬ 
ined positional offsets and uncertainties of a subset of the bright¬ 
est Planck-N\.A sources (38 sources with 70 GHz flux density 
>1.5 Jy, 35 sources with 30 GHz flux density > 1.9 Jy). As ex¬ 
pected, the positional uncertainty was reduced slightly (e.g. from 
59" to 42" at 70 GHz). The offsets in equatorial and ecliptic co¬ 
ordinates changed by 0.5 to 1.0 cr for these bright sources. In 
particular the large offset in Galactic longitude at 30 GHz is re¬ 
duced from 16" to 7" and the large offset in ecliptic latitude at 
100 GHz is reduced from -14" to -7". 


3.3.2. Internal tests 

The positional accuracy of individual sources in the PCCS2 de¬ 
pends on the accuracy of the positions in the detection pipelines. 
This may be evaluated by the injection of sources into the real 
maps. The recovered positions are compared against the known 
positions of the injected sources. Table 6 shows the resulting es¬ 
timates of positional accuracy. Table 7 shows the same thing, 
but here the population of injected sources was limited to those 
with S/N > 20. Note that all of the positional errors are less than 
the width of one HFI map pixel, which is half the width of an 
LFI map pixel. Also note that there is good agreement between 
the positional errors found from the VLA study and those found 
from these simulations. For the lower frequencies the position er¬ 
rors for sources below the 100 % completeness limit must be in¬ 
cluded in the average to bring the position errors for these chan¬ 
nels into agreement with those found from the VLA study. 
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Table 5. Positional offset between the PCCS2 and the VLA (in the sense PCCS2 - VLA) in equatorial, Galactic, and ecliptic 
coordinates. The mean offsets in longitude and latitude (or RA and Dec.) are given, together with the standard error on the mean. 
These offsets place strong limits on any global positional offset between Planck and the VLA, and hence the international celestial 
reference frame (ICRF). The radial positional uncertainty for an individual source, given in the positional error column, is evaluated 
from the standard deviations of the offset in the latitude and longitude positions, assuming no correlations between these offsets. 


Channel 

Equatorial offsets 

Galactic offsets 

Ecliptic offsets 

Positional 

error 

[arcsec] 

Number 
of sources 

RA 

[arcsec] 

Dec 

[arcsec] 

longitude 

[arcsec] 

latitude 

[arcsec] 

longitude 

[arcsec] 

latitude 

[arcsec] 

30 . . . 

-9 + 6 

-7 + 7 

16 + 7 

0 + 6 

-1+7 

-3 + 8 

78.0 

66 

44 . . . 

-11 + 8 

8 + 9 

18 + 9 

-7 + 8 



97.6 

67 

70 . . . 

-3 + 5 

-1 + 6 

4 + 5 

-7 + 6 

-2 + 5 

6 + 6 

58.9 

70 

100 .. . 

-1 + 5 

-4 + 5 

7 + 5 

7 + 5 

9 + 4 

-14 + 5 

53.7 

70 

143 .. . 

-1+5 

-9 + 4 

-6 + 4 

5 + 5 

-1+5 

-13 + 4 

48.8 

66 

217 .. . 

5 + 4 

-1+4 

-2 + 4 

3 + 4 

3 + 4 

-1 + 5 

37.5 

48 


Additionally, we fitted the following functional form relating 
the position error, to the S/N of the detection: 


(rl = 


FWHM\^ 


(7) 


where the value used for the FWHM was the fitted FWHM from 
Table 2, while the values of the parameters ctq and d are shown 
in Table 8. 


3.4. Photometry 

The photometric accuracy of Planck is very high (see Planck 
Collaboration I 2016; Planck Collaboration II 2016; Planck 
Collaboration VII 2016). The consistency of the Planck calibra¬ 
tion is shown in Planck Collaboration I (2016) to be approxi¬ 
mately 0.2 % in most Planck bands. The calibration of Planck, 
however, is based on the measurements of a dipole signal, and 
it is appropriate to ask if that accuracy extends to much smaller 
angular scales. The calibration depends on our knowledge of the 
instruments and the window functions, which in turn depends 
on our understanding of the beam properties in each Planck 
band. Both beam properties and calibration can also be tested by 
comparing, on a statistical basis, the flux densities of compact 
sources at different Planck frequencies or using different photo¬ 
metric methods. We refer to these tests, among others, as “in¬ 
ternal.” We have also undertaken a direct comparison of PCCS2 
flux densities with ground-based or other observations of bright 
sources. We refer to such comparisons as “external” tests. 


3.4.1. Internal Consistency 

Simulations For the HFI channels we characterize the accuracy 
of source photometry by comparing the native flux-density esti¬ 
mates (DETFLUX) of matched sources to the known flux densi¬ 
ties of sources injected into the real maps. The photometric accu¬ 
racy is a function of S/N, and the faint detections are affected by 
upward bias due to noise fluctuations. In the previous PCCS, at 
the higher HFI frequencies, the DETFLUX estimates were found 
to be biased low (Planck Collaboration XXVIII 2014). This has 
been corrected in the construction of the PCCS2. Tables 6 and 7 
show the DETFLUX bias per channel as well as the standard de¬ 
viation of As/cTs, the difference between the input and recovered 
flux densities normalized by the uncertainty on the flux density, 
which would be unity for Gaussian noise. 


Table 6. Native photometry (DETFLUX) bias, (A^), photomet¬ 
ric recovery uncertainty, and radial position uncertainty. The ra¬ 
dial position uncertainty is the 63 % error radius. All were de¬ 
termined from source injection into the maps, using only those 
injected sources with input flux above the 100 % completeness 
threshold. In the three lowest channels, the DETFLUX photom¬ 
etry provided in the catalogues has been corrected for this bias. 


Channel 

DETFLUX bias 
[%] 

stdevCA^/cr^) 

Position error 
[arcsec] 

30 ... . 

-2.34 

0.33 

50.49 

44 ... . 

-4.12 

1.67 

59.57 

70 ... . 

-12.05 

3.69 

44.07 

100 ... . 

1.10 

1.22 

51.96 

143 ... . 

-0.91 

1.44 

43.68 

217 ... . 

-2.36 

1.82 

39.94 

353 ... . 

-3.72 

1.85 

39.59 

545 ... . 

-1.59 

2.13 

39.58 

857 ... . 

-3.51 

2.51 

39.41 


Table 7. As Table 6, but for all detections with S /N > 20. 


Channel 

DETFLUX bias 
[%] 

stdev(As/cr^) 

Position error 
[arcsec] 

30 ... . 

-2.35 

1.16 

37.88 

44 ... . 

-3.15 

1.98 

44.35 

70 ... . 

-13.75 

11.39 

39.69 

100 ... . 

0.58 

1.45 

45.80 

143 ... . 

-1.18 

1.76 

39.53 

217 ... . 

-2.06 

2.15 

38.33 

353 ... . 

-3.24 

2.14 

38.57 

545 ... . 

-0.81 

2.54 

37.85 

857 ... . 

-2.27 

2.79 

37.99 


Comparisons of the four different flux-density estimates We 
next compare values derived from the four different methods 
of assessing flux densities. Figures 10, 11, 12, and 13 show 
the results for four Planck channels. These comparisons are 
made against APERFLUX. This flux estimation method is the 
simplest and makes the smallest number of assumptions about 
the data. However, DETFLUX has smaller uncertainties than 
APERFLUX, which may be seen in these figures by the up¬ 
ward curve towards lower APERFLUX values in the compar¬ 
ison with DETFLUX. This may be understood as there being 
a clear signal present in DETFLUX when the APERFLUX is 
compatible with noise. Figure 14 shows the comparison against 
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Table 8. Parameters (Tq and d determined by fitting Eq. (7), re¬ 
lating the position error, cr^, from the simulations to the S/N of 
the detection. 


Channel 

o-Q 

[arcmin] 

d 

30 ... . 

0.267 + 0.001 

2.14 + 0.44 

44 ... . 

0.217 + 0.001 

1.74 + 0.59 

70 ... . 

0.538 + 0.001 

1.68 + 0.23 

100 ... . 

0.685 + 0.001 

1.61 + 0.01 

143 ... . 

0.615 + 0.001 

1.60 + 0.01 

217 ... . 

0.580 + 0.002 

1.38 + 0.02 

353 ... . 

0.578 + 0.002 

1.43 + 0.03 

545 ... . 

0.539 + 0.002 

1.48 + 0.04 

857 ... . 

0.546 + 0.0004 

1.46 + 0.02 



APERFLUX [Jy] 


Fig. 10. Comparison of the DETFLUX, PSFFLUX, and 
GAUFLUX fiux-density estimates with APERFLUX for the 
PCCS2 30 GHz catalogue. The fractional difference is defined 
as {S - 5 aperflux)/^'APERFLUX. The green and blue points cor¬ 
respond to sources where 5'aperflux/‘5'aperflux_err > 5. Grey 
and green points correspond to sources with \h\ < 5°, while the 
red and blue have |Z?| > 5°. 


DETFLUX for APERFLUX at 30, 70, and 143 GHz. In these 
plots, as expected, the curvature disappears, and we see good 
agreement between the methods that becomes progressively 
noisier towards lower values of DETFLUX. For unresolved 
sources in regions where there is little non-Gaussianity present 
in the background, DETFLUX is the flux estimation method 
of choice, given its greater sensitivity. However, in regions of 
high non-Gaussian background emission, DETFLUX is less ro¬ 
bust. This may been seen by the lack of consistency between 
DETFLUX and APERFLUX, in these figures, at least for the 
green and grey points representing sources that lie within 5° 
of the Galactic plane. As the frequency increases so do the 
levels of non-Gaussian emission. In the comparison between 
DETFLUX and APERFLUX at 353 GHz, which only contains 
sources at the higher Galactic latitudes, there is a large degree 
of scatter. Indeed, at 353 GHz and above it is advisable to favour 
APERFLUX over DETFLUX. 



APERFLUX [Jy] 


Fig. 11. Comparison of the DETFLUX, PSFFLUX, and 
GAUFLUX fiux-density estimates with APERFLUX for the 
PCCS2 70 GHz catalogue. The fractional difference is defined 
as {S - 5 aperflux)/5'APERFLUX- The green and blue points cor¬ 
respond to sources where 5'aperflux/5'aperflux_err > 5. Grey 
and green points correspond to sources with \b\ < 5°, while the 
red and blue have |Z?| > 5°. 



APERFLUX [Jy] 


Fig. 12. Comparison of the DETFLUX, PSFFLUX, and 
GAUFLUX fiux-density estimates with APERFLUX for the 
PCCS2 143 GHz catalogue. The fractional difference is defined 
as {S - 5*APERFLUX)/5'APERFLUX- The blue points correspond to 
sources where 5'aperflux/5'aperflux_err > 5. 


Interband Comparisons As an additional internal test, we com¬ 
pared PCCS2 flux densities at one band with those in neigh¬ 
bouring bands. We performed this comparison for the six low¬ 
est Planck channels. We began by selecting all PCCS2 sources 
at 70 GHz with the following restrictions: flux density 5(70) > 
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APERFLUX [Jy] 


Fig. 13. Comparison of the DETFLUX, PSFFLUX, and 
GAUFFUX flux-density estimates with APERFFUX for the 
PCCS2 353 GHz catalogue. The fractional difference is defined 
as (S - 5aperflux)/5'APERFLUX- The blue points correspond to 
sources where 5'aperflux/5'aperflux_err > 5. 



DETFLUX [Jy] 


Fig. 14. Comparison of the APERFFUX flux-density esti¬ 
mates with DETFFUX for the PCCS2 30 GHz (top panel), 
70 GHz (middle panel) and 143 GHz (bottom panel) cata¬ 
logues. The fractional difference is defined as (^aperflux “ 
5'DETFLUx)/^'DETFLUX- Grey point are sources at \b\ < 5°; while 
red points have \b\ > 5°. The greater depth of the 143 GHz cata¬ 
logue is clear. 


0.9 Jy (virtually all such sources had S/N >7); Galactic latitude 
\b\ >10°; and no evidence of extension. We then matched these 
with sources in the PCCS2 at 30, 44, 100, 143, and 217 GHz. 
Of the 203 sources, more than 99 % were detected at 30, 44, and 
100 GHz as well as 70 GHz, and 97% at 143 GHz. Since vir¬ 
tually all the sources had synchrotron spectra, and, in general. 


Table 9. Slope found in the plots of predicted versus measured 
flux densities in five of the Planck bands. A slope > 1 im¬ 
plies that the measured flux density exceeds the predicted value. 
Departures from unit slope can be explained by slight curvature 
of the spectra; we show the required change in spectral index, 
Aa, to explain the departure from unity. 


Channel 

Slope 

Required Acr 

44 ... . 

0.985 + 0.005 

+0.066 

70 ... . 

1.025 + 0.005 

-0.125 

100 ... . 

1.038 + 0.005 

-0.210 

143 ... . 

0.997 + 0.005 

+0.015 

217 ... . 

1.017 + 0.005 

-0.076 


negative spectral indices, it is not surprising that only 90 % of 
the sources could be identified at 217 GHz. We note that some 
of the sources not found in the PCCS2 at 217 GHz did appear in 
the PCCS2E; these few sources, however, were not used in this 
internal test. 

We used this merged catalogue to calculate spectral indices 
for each of the 203 sources. The spectral indices were then used 
to make the small colour-corrections to the flux density in each 
band for each source. For the colour-corrections at 30 GHz we 
used the 30-44 spectral index, and for 353 GHz the 217-353 
spectral index. For the other five bands, for band A, we used 
the spectral index found between bands N - \ and N + 1 (e.g. 
for 44 GHz, we used the 30-70 spectral index). The colour- 
corrections we used are tabulated in Planck Collaboration II 
(2016) and Planck Collaboration VII (2016).The amplitude of 
these corrections in a given Planck band ranges from around 
2.5 % at 70 GHz to less than 1 % at 30 GHz. For a given source, 
the precision of the colour-corrections was typically 0.2 %. 

Next, we used the colour-corrected flux densities to recom¬ 
pute spectral indices for each source. These spectral indices were 
then used to predict a flux density for each source in frequency 
band N by assuming a constant spectral index between bands 
N - \ and N + 1. For instance, we interpolated between the 
28.4 and 70.4 GHz flux densities to predict a 44.1 GHz flux den¬ 
sity for each source, using the calculated spectral index for that 
source. These predicted values were then plotted against the ac¬ 
tual (colour-corrected) measurements at 44.1 GHz. This opera¬ 
tion was repeated for 70, 100, 143, and 217 GHz. If the flux- 
density scales of Planck are consistent across bands and if the 
spectral index is constant as assumed, we expect to see lines of 
unit slope. In fact, the slopes were close to unity for all five bands 
tested, as shown in Table 9. 

We examined several different possibilities for the slight de¬ 
partures from unit slopes. First, we explored the possibility that 
CO line emission could influence the results by perturbing the 
flux densities. The PCCS2 flux densities used were the MHW2 
estimates. The Altering removes all foreground emissions on 
large scales; only scales approaching the size of the beam could 
affect the MHW2 flux density estimates. Thus, Galactic CO 
emission could introduce scatter in the 100 GHz values, but 
should not give significant coherent offsets. There is, however, 
the question of redshifted CO line emission from the sources 
themselves. Since the sources meeting our selection criteria are 
mostly bright blazars (Planck Collaboration XIV 2011), the ra¬ 
tio of CO line flux to continuum emission is expected to be very 
small. 

Next, we investigated whether the statistically significant de¬ 
partures from unit slope in the plots of predicted versus mea¬ 
sured flux could reasonably be explained by a breakdown in the 
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Table 10. Comparison of the corrected Planck flux-density mea¬ 
surements with the ATCA, VLA, and ACT values. The error in 
the flux-density ratios is purely due to statistical uncertainties; 
the error in the % increment includes the uncertainties due to the 
Planck beams and calibration. 


Frequency 

[GHz] 

Flux density ratio 
(ground/P/afzcG 

% increment of Planck over 
ground-based measures 

22.45 (ATCA) . . 

0.99 + 0.017 

1 + 1.7 

28.45 (VLA) . . . 

0.97 + 0.008 

3 + 1.2 

43.34 (VLA, ATCA) 

0.94 + 0.013 

6 + 1.7 

147.70 (ACT) . . . 

0.97 + 0.03 

3 + 3 

217.60 (ACT) . . . 

0.96 + 0.03 

4 + 3 


assumption that the spectral index of the sources stays constant 
from band A - 1 to band N +1. For instance, the measured slope 
at 70 GHz is 1.0125; that is, the measured 70 GHz flux densi¬ 
ties are about 2.55 % higher than we would And by interpolating 
between 44 and 100 GHz, assuming no change in the spectral 
index between 44 and 100 GHz. If instead we allow for spectral 
curvature, the small discrepancy can be reduced or eliminated. 
With the simplest assumption, a sharp change of spectral index, 
Aa, at 70 GHz, we And that the small excess of measured over 
predicted flux can be explained by Aa = -0.125 at 70 GHz. We 
performed similar calculations for this simple model in the other 
Planck bands (see Table 9). Spectral index changes of this mag¬ 
nitude are reasonable (electron ageing can account for a change 
in spectral index of around -0.5); see Sect. 4 for plots of spectral 
index distributions for all sources in the PCCS2, not just those 
used in this analysis. 

3.4.2. External consistency 

The calibration of Planck is precise, and we have demonstrated 
the internal consistency of the flux densities in Sect. 3.4.1. The 
calibration of Planck is also absolute in the sense that it de¬ 
pends only on the motion of the satellite and the 0.02% ac¬ 
curate measurement of the CMB temperature (Fixsen 2009). 
Consequently, comparing Planck flux densities to those mea¬ 
sured by other instruments is actually a check on the accuracy 
of the latter. Indeed, Butler et al. (2015) have employed Planck 
measurements to reflne the centimetre-wavelength flux-density 
scales used at the VLA and the ATCA. Here, we summarize the 
results of that study and of comparisons with other CMB and 
submillimetre instruments and missions. 

The 30 and 44 GHz channels The comparison between Planck 
flux densities at 28.4 and 44.1 GHz and ground based observa¬ 
tions at 22.45, 28.45, and 43.34 GHz is based on observations 
carried out at the ATCA and the VLA in April and May 2013. 
Both instruments observed a set of strong, unresolved, uncon¬ 
fused radio sources also scanned by Planck in this time inter¬ 
val. These observations are part of a wider effort by Perley and 
Stevens (2015, in preparation) to compare the flux-density scales 
of the two interferometers, the VLA in the north and ATCA in 
the south. The Planck DETFLUX measurements were derived 
not from the PCCS2, which averages over the four years of LFI 
observations, but from a special map constructed using only data 
from 1 April to 30 June 2013. This was necessary in order to 
minimize the effects of source variability. Since the central fre¬ 
quencies of the Planck bands did not exactly match the frequen- 




Planck flux density [Jy] 

Fig. 15. Comparison of simultaneous colour-corrected flux den¬ 
sity measurements by Planck and the VLA at 28.45 GHz (top) 
and Planck and the VLA (dots) or ATCA (open squares) 
43.34 GHz (bottom). The best fit line is shown in red. The Planck 
44 GHz channel is noisier than the 30 GHz channel and shows 
Eddington bias at low flux densities. Some sources have larger 
error bars than others because they lie close to the Galactic plane 
where the uncertainties in the measured flux densities can be 
large. The tiny VLA error bars have not been plotted. 


cies employed by the ground-based instruments, we interpolated 
and colour-corrected the Planck measurements to 22.45, 28.45, 
and 43.34 GHz. For both purposes, we used spectral indices de¬ 
rived from the far more precise interferometric measurements 
(VLA or ATCA). A comparison of the corrected Planck mea¬ 
surements with the ATCA and VLA values is shown in Table 10 
and Fig. 15. The Planck measurements are consistently slightly 
higher across all frequencies. These results are summarized in 
Butler et al. (2015), and described in greater detail in Partridge 
et al. (2015) where various tests of the validity of the results are 
presented. In particular, the effect of dropping sources found to 
have varied over the three-month period of the Planck observa¬ 
tions is examined. The results in Table 10 are estimates derived 
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from Butler et al. (2015). The discrepancies between the satel¬ 
lite and ground-based values lie close to or within the estimated 
error in the latter. For the flux-density scale of Perley & Butler 
(2013), employed at the VLA, this uncertainty is estimated to be 
5 %, and roughly the same level of precision may be assigned to 
the flux-density scale employed at ATCA. We also compared the 
ATCA and VLA measurements to PCCS2 30 and 44 GHz flux 
densities (averaged over 4 years) and found consistent results, 
for the same sources. As expected, source variability substan¬ 
tially increased the scatter. Excluding three manifestly variable 
sources, we And Planck = 0.96 + 0.02 at 28.45 GHz and 
{Nl^AScKYCAWlanck = 0.93 + 0.03 at 43.34 GHz. These are 
consistent with the more precise comparison described above. 
The 3% and 6% differences may be compared to the quoted 5% 
uncertainty in the flux density scales as given by Perley & Butler 
(2013). 

The Metsahovi Observatory is continuously monitoring 
bright radio sources in the northern sky at 37 GHz (Terasranta 
et al. 2004). From their sample, sources brighter than 1 Jy were 
selected and their flux densities averaged over the period of 
Planck observations used for the PCCS2 (Planck Intermediate 
results, in preparation); note that this period corresponds to the 
full duration of the Planck mission. Hence, the uncertainties in 
Fig. 16 reflect the variability of the sources during the Planck 
mission. The Planck measurements were colour-corrected and 
extrapolated to the Metsahovi frequency before making the com¬ 
parison. The Planck and Metsahovi flux densities agree at the 
0.3 % and 0.1% level, and an uncertainty of +4 %, at 30-44 GHz 
and 30-70 GHz, respectively. 

The 143 and 217 GHz channels Flux densities at 143 and 
217 GHz from the earlier PCCS were compared to measure¬ 
ments made at the Atacama Cosmology Telescope (ACT) by 
Louis et al. (2014). We repeated those comparisons using the 
new PCCS2 values, and extended the comparison to South Pole 
Telescope (SPT) flux densities from Mocanu et al. (2013b). We 
employed the DETFLUX values from the PCCS2. As in Louis 
et al. (2014), the Planck flux densities were colour-corrected and 
extrapolated to the central frequencies of both ACT (147.6 and 
217.6 GHz) and SPT (152.9 and 218.1 GHz), using the spectral 
index appropriate for each source. It is important to bear in mind 
that the ground-based measurements were in almost all cases 
far from simultaneous with those of Planck. Thus the variabil¬ 
ity of sources, virtually all AGN, caused significant scatter. In 
the case of the ACT equatorial sources, we had measurements 
from both the 2009 and the 2010 seasons. That allowed us to And 
and drop two manifestly variable sources. We also dropped two 
sources with thermal spectra. For SPT, Mocanu et al. (2013b) 
present just a single flux density for each source, so in most 
cases we had no means of discovering and removing sources 
that were variable. Four sources, however, were detected by all 
three experiments, ACT and SPT in 2008 and Planck integrated 
over the mission. Three of these were evidently variable (and 
are discussed further below). For the comparison between cor¬ 
rected Planck flux densities and those of ACT at 147.6 GHz, 
we were left with 58 sources in common. ACT flux densities 
were on average 0.97+0.03 times Planck's. At 217.6 GHz, fewer 
ACT sources (50) were detected by Planck, and we And ACT = 
0.89+0.03 times Planck. 

We now consider the effect on these results of sources known 
to be variable. First, comparison of ACT measurements made 
in 2009 and 2010 showed that two sources varied strongly. 
Dropping them changed the 147.6 GHz slope to 0.95. On the 




Metsahovi flux density [Jy] 


Fig. 16. Comparison between the Metsahovi and the colour- 
corrected PCCS2 flux densities (DETFLUX) interpolated to 
37 GHz using 30 and 44 GHz {top) and 30 and 70 GHz {bottom). 
The multiple observations of each source have been averaged to 
a single flux density; the averaging was performed over the pe¬ 
riod of the full Planck mission, not just the epochs at which each 
source was observed by Planck. The uncertainties, therefore, re¬ 
flect the variability of the sources instead of the flux-density ac¬ 
curacy of the measurements, which is of the order of a few mJy. 


other hand, if we drop the three variable sources detected by 
all three experiments, the slope becomes 0.98. If we drop all 
five sources for which we have direct evidence of variability, 
the slope at 147.6 GHz settles to 0.97+0.02. Since the exclusion 
of variable sources moves the slope both up and down in am¬ 
plitude, we adopt 0.97+0.03 for the relation between ACT and 
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Fig. 17. Comparison of flux densities from Planck and Herschel 
at 350yum (857 GHz) using the Herschel Reference Survey. The 
one-to-one lines are shown in red. 


Planck flux densities; ACT flux densities are about 3 % lower 
than Planck's. At 217.6 GHz, if we remove the few sources for 
which we have direct evidence of variability, the slope changes 
to 0.96+0.03. 

Planck found fewer SPT sources than ACT sources (25 at 
152.9 GHz and 30 at 218.1 GHz), and the scatter due to source 
variability was larger. If we include all sources, we again find the 
ground-based fiux densities are lower than Planck's: 0.95+0.05 
at 152.9 GHz and 0.88+0.05 at 218.1 GHz. If we now exclude 
the three sources that were seen to vary between 2008 and 
the later Planck mission, the SPT results become 0.99+0.05 at 
152.9 GHz and 0.97+0.05 at 218.1 GHz. We adopt these values, 
all lying between 0.96 and 0.99, as evidence of good agreement 
between the ground-based and Planck fiux density scales at 143 
and 217 GHz. The small differences between the Planck fiux 
density scales at these frequencies and those measured from the 
ground can be compared to the following uncertainties quoted 
for the ground based experiments: SPT, 1.6% and 2.4% calibra¬ 
tion uncertainty in the maps at 150 and 220 GHz, respectively 
(Mocanu et al. 2013b); and ACT, 6% uncertainty in the fiux den¬ 
sity scale at 148 GHz (Louis et al. 2014). 

The 857 GHz channel At this frequency, there are no other all¬ 
sky surveys available matching in frequency against which the 
PCCS2 fluxes can be compared. However, the 350yum chan¬ 
nel of the SPIRE instrument (Griffin et al. 2010) on Herschel 
(Pilbratt et al. 2010) is a close match to the 857 GHz passband 
of Planck. Pointed observations of compact extragalactic tar¬ 
gets obtained with SPIRE can thus be used to validate the fiux 
density measurements made in the PCCS2. The most useful set 


of observations from SPIRE for our purposes are those from 
the Herschel Reference Survey (HRS) (Boselli et al. 2010), a 
survey of the far-infared and submillimetre properties of local 
bright galaxies. The published SPIRE photometry for this sam¬ 
ple (Ciesla et al. 2012) uses either PSE fitting for sources unre¬ 
solved by Herschel or apertures matched to the observed sizes of 
the sources in the Herschel maps. In both these cases the aper¬ 
tures used will be much smaller than the Planck beam. Since 
several of the HRS galaxies have nearby bright companions, and 
since other sources might be included in the large Planck beam, 
new fiux density values were extracted for the HRS sources us¬ 
ing apertures matched to the size of the Planck beams (Eales 
et al., private communication). This allows a direct compari¬ 
son of the SPIRE 350yum fiux densities to the Planck 857 GHz 
fluxes for the same 141 objects. It should be noted that only 
three of these objects have not been flagged as EXTENDED by 
Planck. We should, therefore, not be surprised if the fiux ex¬ 
traction methods (that assume a single point-like source) will be 
biased low. In Eig. 17 and Table 11 we show the results of this 
comparison for the four different flux extraction methods used 
in the PCCS2; namely DETELUX, APERELUX, PSEELUX, 
and GAUELUX. The best performing Planck flux-density ex¬ 
traction method is, perhaps not surprisingly, the method that 
most closely resembles the flux-density extraction method ap¬ 
plied to the Herschel maps, namely APERELUX, which shows 
good agreement between Herschel and Planck fluxes over the 
full range of source brightness. The worst performing method, 
in contrast, is DETELUX, which shows an increased scatter, and 
an overall bias to lower Planck values for the brighter sources. 
This bias is expected if the sources are not truly point-like. Eor 
a population of slightly extended sources we expect a noticeable 
downward bias for the brighter sources, which disappears into 
the noise for fainter objects. This pattern of bias for the brighter 
sources which disappears once these sources are excluded is also 
seen for PSEELUX. The increased scatter seen for the brighter 
sources with DETELUX probably arises from errors in the re¬ 
covered position. This is due to the relationship between the 
scale of the wavelet used and the pixel size in the maps. Eor 
this channel (857 GHz) the small beam size and the increased 
level of foregrounds in the map mean that the optimum scale 
for the wavelet is very narrow with respect to the map pixels. 
This makes the flux-density estimate extremely sensitive to er¬ 
rors in the recovered position of the source. This effect is most 
pronounced at this channel and for the brighter sources. A sim¬ 
ilar bias seen in GAUELUX is a little more difficult to explain, 
but again may result from the assumption that sources are single 
and point-like. A double source, for instance, may pose diffi¬ 
culties in the flux-density estimation resulting in a bias low. We 
conclude that there is a good match between the PCCS2 fiux 
densities at 857 GHz and Herschel fiux densities in the match¬ 
ing 350/im SPIRE band. We also conclude that, for most pur¬ 
poses concerned with the flux-density measurement of compact 
sources like the galaxies discussed here, APERELUX is the most 
appropriate flux-density measure to use for the higher frequency 
channels. 

3.5. Polarization measures 

To validate the polarization measurements in the PCCS2 and 
PCCS2E catalogues, we rely on simulations based on the injec¬ 
tion of sources with known properties into the Planck polariza¬ 
tion maps between 30 and 70 GHz and the EEP8 simulated maps 
(Planck Collaboration XII 2016) between 100 and 353 GHz. It 
should be noted that the injection of sources into the real maps 
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Table 11. Comparison of flux-density measurements from 
Planck and Herschel at 857 GHz (350 jam). The Planck measure¬ 
ments for all flux-extraction methods are lower than the Herschel 
results. If, however, the 17 sources in the HRS sample with flux 
densities greater than 10 Jy are excluded from the analysis, then 
the flux-density ratios for all flux-extraction methods are within 
1 (T of unity. 


Flux density ratio {Herschel/Planck) 


method All 141 HRS sources Discarding S'hrs >10 Jy 


DETFLUX. 1.045 + 0.010 0.982 + 0.019 

APERFLUX_ 1.019 + 0.011 1.009 + 0.021 

PSFFLUX. 1.065 + 0.022 1.016 + 0.038 

GAUFLUX. 1.056 + 0.013 1.024 + 0.024 


has the advantage that many thousands of sources can be used 
to test the analysis; however this procedure does not simulate 
the leakage due to the bandpass mismatch. The FFP8 maps do 
contain this effect, although there are only hundreds rather than 
thousands of polarized sources. Both approaches, however, al¬ 
low us to test the fidelity of the polarized flux densities and po¬ 
larization angles produced by our analysis pipelines. They also 
allow us to compute measures for the completeness and reliabil¬ 
ity of recovered polarization measurements. Since the extraction 
of these measurements is non-blind, based on the positions pro¬ 
vided by the analysis of the temperature maps, these terms for 
the polarization pipelines only have any meaning given that the 
source is real. 

We define the “polarization photometric completeness” 
(PPC) as the percentage of polarized sources correctly identified 
as polarized above a given true polarized flux density, and “po¬ 
larization photometric reliability” (PPR) as the percentage of po¬ 
larized sources whose polarized flux density is contained in the 
interval defined by the best-fit value and +3 cr errors. Figure 18 
shows the PPC and PPR as a function of the true and estimated 
polarized flux density, respectively, for the PCCS2 for all the 
polarized channels. The red dashed lines show the results for 
the subsets of significantly polarized sources, whose measure¬ 
ments are provided for both LFI and HFI, while the green dot- 
dashed lines show the results when marginally polarized sources 
are also included in the analysis. Here we see that the inclusion 
of the marginal data increases the completeness with only a min¬ 
imal decrease in the 353 GHz channel reliability. The negative 
“kinks” in the reliability curves of the 217 and 353 GHz channels 
are caused by the same single source. This source was detected 
in the intensity maps with S/N =151.5 (217 GHz) and the recov¬ 
ered flux density was underestimated by around 23 cr, where cr 
is the estimated error on the flux density. The dramatic underes¬ 
timation of the flux density was caused by the recovered position 
of the source being offset from the true position by 1.27 arcmin. 
In this S/N regime such an offset is sufficient to explain the fail¬ 
ure of the PPR criterion, as the polarized flux density will be 
underestimated by more than 3 cr. 

A comparison between the LFI and HFI implementations 
of the common procedure for extracting the polarization mea¬ 
surements for the significantly polarized sources was performed. 
Figure 19 shows a comparison of the recovered polarized flux 
densities (for each implementation), and the true polarized flux 
density. We have assessed the performance of each method, and 
And that above 250 mJy, where we are complete, the average re¬ 
covered polarized flux density is within 1 % of the true value for 
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Fig. 18. Left: polarized photometric ompleteness (PPC). Right: 
polarized photometric reliability (PPR). Red dashed lines: de¬ 
rived with the common method applied to both LFI and HFI. 
Green dot-dashed lines: including the marginal polarization data 
(HFI only). From top to bottom, 30-353 GHz. The plots were 
constructed using the PCCS2catalogue only. 


IFCAPOL and within 0.8 % for PwSPOL, and that they are within 
0.2 % of each other. 


3.5.1. Internal consistency 

30-70 GHz In order to assess the performance of the maxi¬ 
mum likelihood estimator used to provide the polarization mea¬ 
surements in the three lowest Planck channels (IFCAPOL), we 
used Monte Carlo simulations. Point sources were simulated 
and convolved with the appropriate Gaussian effective beam 
from Table 2 for each channel, and were injected into HEALPix 
A^side = 4096 maps. These maps were then downgraded to 
A^g-^g = 1024 to match the pixelization of the LFI maps, and 
the source maps were added to the Planck 2015 (PR2) Q and U 
maps. The sources were injected away from known bright radio 
sources. In all, 37 000 sources, in 50 flux-density bins, were in¬ 
jected into the Q and U maps. The pipelines were run given the 
positions of these simulated sources, producing polarized flux 
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Fig. 19. Comparison between the LFI and HFI implementations 
of the maximum-likelihood estimator code to measure polar¬ 
ized flux densities. The comparison was performed using the 
100 GHz FFP8 maps, and the locations of the point sources in¬ 
cluded in the simulation, because the polarization pipelines per¬ 
form non-blind extractions. The recovered fractional difference 
of the polarized flux densities of the two methods, defined as 
(recovered - true)/true, is plotted against the true polarized flux 
density, as simulated in the FFP8 maps. The results from the LFI 
pipeline, IFCAPOL, are shown by the green diamonds, and the 
results from the HFI pipeline, PwSPOL, are shown by the blue 
crosses. 


density and polarization angle measurements. The results of the 
Monte Carlo simulations on the recovery of the polarized flux 
densities are shown in Fig. 20. We can see that the polarized flux 
densities are recovered in an unbiased way for strongly polarized 
sources, and that the faintest ones suffer from Eddington-type 
bias. Figure 21 shows the equivalent plot for the polarization 
angle and here we can see that the angle is recovered success¬ 
fully for all channels, where the uncertainties increase towards 
fainter polarized flux densities. Based on these simulations, we 
computed the polarization photometric completeness and relia¬ 
bility, shown in Fig. 18. In the simulations we injected sources 
at all Galactic latitudes, and we did not apply any Galactic cut. 
We And that our catalogues are complete at the 90 % level at 
200 mJy (polarized flux density) at 30 GHz, and at 400 mJy at 44 
and 70 GHz. At the 600-700 mJy level the three catalogues are 
complete. A summary of these results can be found in Table 14 
(Sect. 4). 


100-353 GHz The recovery of the polarized flux densities from 
the FFP8 simulations, using the method common to both LFI 
and HFI, is shown in Fig. 22 for the PCCS2 and the PCCS2E. 
The PCCS2 polarized flux density estimates are reliable across 
the full range of flux densities and considerably more reliable 
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Fig. 20. Internal validation using Monte Carlo simulations to as¬ 
sess the recovery of the polarized flux density for the 30, 44, 
and 70 GHz channels. In these simulations, point sources were 
injected into the Q and U maps with 50 different polarized 
flux density values, starting at 0.2 Jy and increasing with a step 
size of 0.1 Jy. The fractional difference, defined as (recovered - 
true)/true, is plotted against the true polarized flux density. The 
recovered polarized flux densities are unbiased, except for the 
faintest sources, were the effect of Eddington-type bias is seen. 


than the PCCS2E. Given that the PCCS2E catalogue covers the 
Galactic plane region this is hardly unexpected. The polarized 
flux densities are unbiased over the range of values where the 
survey is complete, but for the fainter sources a positive bias is 
present, as expected from Eddington-type bias. Figure 23 shows 
the recovery of the polarization angle, which is unbiased over 
the full polarized flux-density range, although the errors in its 
recovery increase as the polarized flux density of the source de¬ 
creases. This behaviour is exactly what is expected if our as¬ 
sumptions made in Sect. C.2 hold; namely that (Tq ^ cru, and 
that they are uncorrelated. Consider sources with a given polar¬ 
ization angle in the range of polarized flux densities which suf¬ 
fer from Eddington bias: the required upward noise fluctuations 
could come from 2 or f/ or both. Once the polarization angle is 
evaluated, if our assumptions hold, then the average value found 
for the angle will be unbiased but its dispersion will be much 
larger than for the brighter sources. In Fig. 23 we also see that 
the measurements in the PCCS2E are again less reliable than the 
PCCS2. 


3.5.2. External consistency 

The limited number of polarimetric millimetre surveys, the small 
number of bright Planck sources with a high significance in po¬ 
larization, and the fact that the majority of polarized sources are 
variable, makes it difficult to validate the Planck polarized flux 
densities with external datasets. 

One of the objects from the catalogue that we have studied 
in detail is Tau A, also known as Ml or the Crab Nebula. This 
object is resolved in the higher frequency Planck channels and 
may not be the best source for validation at these frequencies; 
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Fig. 21. Internal validation using Monte Carlo simulations to as¬ 
sess the recovery of the polarization angle of the 30, 44, and 
70 GHz channels. In these simulations, point sources spanning 
the full range of polarization angles were injected into the Q and 
U maps with 50 different polarized flux density values, starting 
at 0.2 Jy and increasing with a step size of 0.1 Jy. The difference 
between the recovered and true angles is plotted against the true 
polarized flux density. Here we see there are no biases in the 
recovery of the polarization angle, although obviously the un¬ 
certainty increases for fainter sources. 


nevertheless it is the brightest compact source in polarization in 
Planck and has been thoroughly studied in other experiments. In 
Table 12 we compare the total intensity and polarized flux den¬ 
sities, polarization fraction, and polarization position angle for 
TauA with measurements from WMAP (Weiland et al. 2011) 
and the IRAM 30m telescope at 89 GHz (Aumont et al. 2010). 
In general, it is assumed that the polarization position angle of 
Tau A is constant across the frequency range of interest, up to at 
least 353 GHz. In Table 12 one can see that the Planck polariza¬ 
tion position angles (as measured by the maximum likelihood 
Altering method and by aperture photometry) are significantly 
different from those of WMAP and IRAM at some frequencies. 
We have investigated these discrepancies and found that there 
are multiple factors in LFI and HFI affecting our measurements 
of the position angles of Tau A that deserve further attention. 

First, at the position of TauA in the Stokes U maps one 
can see a small spurious signal that is affecting the angle mea¬ 
surements. For this object, most of the polarized signal is in the 
Stokes Q map, and U makes very little contribution to the total 
polarized flux density (quadrature sum of Q and U). However, as 
shown in Table 12, when we calculate the polarization position 
angle with either method, the contribution from the spurious sig¬ 
nal has the effect of changing the position angles by up to 5° from 
the -88?2 measured at IRAM (Aumont et al. 2010). This spuri¬ 
ous signal has been introduced by the complex cross-terms in the 
polarized beams, which are normally expected to be very small 
and below the map noise level, but show up here because Tau 
A is so bright. In an attempt to remove this effect, we have pro¬ 
duced maps where we have deconvolved the beam and fitted for 
spurious signals. Results for LFI frequencies, from the ArtDeco 



Polarized flux density [Jy] 


Fig. 22. Fractional difference between the true and recovered po¬ 
larized flux densities, from top to bottom, for the 100-353 GHz 
FFP8 simulations. Left (red): PCCS2. Right (green): PCCS2E. 
Here the fractional difference is defined as recovered minus true 
divided by the true value. These are the significantly polarized 
sources, as found by the common method, and the uncertain¬ 
ties associated with the best-fit estimates are +3 cr error bars. 
Eddington-type bias is seen at lower polarized flux densities. 


map-making pipeline at the LFI DPC, are shown in column 7 of 
Table 12. These new measures of position angle agree with the 
expected value from IRAM and WMAP. However, we caution 
that these new maps are under development and the measure¬ 
ments should be used with caution. There is an ongoing effort, 
not yet completed, to generate a similar set of maps for the HFI 
frequencies, in order to understand whether similar spurious sig¬ 
nals affect the HFI polarization angle measurements. The tests 
that we have carried out with the new LFI maps indicate that 
the polarized flux densities and angles of the other sources in 
the catalogue are mostly unchanged, but some sources with high 
polarized flux density may be marginally affected. 

Second, the angular size of Tau A measured with IRAM at 
89 GHz (Aumont et al. 2010) shows that it will be slightly re¬ 
solved by Planck in the HFI channels. This could have an impact 
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Polarized flux density [Jy] 


Fig. 23. Difference between the recovered and true polarization 
angles, from top to bottom, for the 100-353 GHz FFP8 simula¬ 
tions. Left (red): PCCS2. Right (green): PCCS2E. These are the 
significantly polarized sources, as found by the common method, 
and the uncertainties associated with the best-fit estimates are 
+3 (T error bars. 


on the flux densities and polarization position angles measured 
with both the filtering method and with aperture photometry. In 
particular, the filtering method assumes that the sources are un¬ 
resolved, so the flux densities derived with this method should 
be regarded as lower limits if the source is extended. In the case 
of aperture photometry, the integration radius that we use in HFI 
is too small for a source as large as TauA. We have therefore 
increased the aperture radii to 2 x (FWHM^ -P where 

^TauA = 4' as measured by IRAM, and recomputed the angles 
for the four HFI channels with polarization capabilities. This im¬ 
proves the agreement between the HFI position-angle measure¬ 
ments and those from WMAP and IRAM, except at 353 GHz. 
The results are shown in column 7 of Table 12. 

Third, as shown in Planck Collaboration VI (2016) and 
Planck Collaboration VIII (2016), different methods have been 
used in LFI and HFI for correcting the frequency maps for a 
global leakage signal from total intensity into polarization due 
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Fig.24. Polarization position angles for TauA from WMAP 
(blue squares, Weiland et al. 201 1) and Planck (black dots). The 
IRAM measurement from Aumont et al. (2010) is also shown as 
a solid line, where a +1(T uncertainty is shown with dashed lines. 


to the bandpass mismatch. In practice, this means that the LFI 
frequency maps have not been corrected, whereas the HFI ones 
have. The global bandpass corrections applied to the maps are 
not the same as the local corrections that we apply in this paper 
for point sources, which depend on the spectrum of each source 
(see Appendix C). However, even though our techniques to ex¬ 
tract the flux densities of the sources in the Q and U maps will 
remove most (if not all) of the global correction if the sources 
are point-like, this may not be true for extended sources. To test 
this, we have extracted the polarization position angles from HFI 
maps where the global bandpass correction had not been applied. 
The results are shown in Table 12 and in Fig. 24, where we com¬ 
pare the polarization position angles from WMAP and IRAM 
with the new measurements from Planck. The recovered angles 
between 100 and 217 GHz are very similar to the case where the 
global bandpass mismatch had been applied, except at 353 GHz, 
where the new angle is in much better agreement with the other 
Planck, WMAP, and IRAM measurements. 

These analyses show the complexity of the Planck maps in 
polarization, particularly for bright extended objects like TauA, 
where we had to fit and remove the spurious signal in the LFI, 
increase the integration aperture radii in the HFI, and remove the 
global bandpass correction at 353 GHz, in order to achieve a con¬ 
sistency within 2° in polarization position angle between Planck 
and WMAP or IRAM. Tau A is the only very bright polarized 
source in our maps, and the limited amount of polarization in¬ 
formation available at Planck frequencies for other sources lim¬ 
its our ability to conclude that the polarization angles of the rest 
of the sources in the catalogue are not affected by these issues. 
Therefore, as in the case of Tau A, the polarization position an¬ 
gles should be used with caution. In future releases of the Planck 
products we will revisit this issue. 

In addition, we have cross-matched the Plateau de Bure in¬ 
terferometer (PdBI) polarimetric survey of 86 AGN at 100 GHz 
(Trippe et al. 2010) with the Planck 70 and 100 GHz PCCS2 
catalogues, finding two sources in common (PKS 0851-r202 and 
3C273). For these sources we see good agreement: the Planck 
polarized flux densities are 509 + 106mJy and 515 + 102mJy 
at 70 GHz, and 566 + 38mJy and 503 + 36mJy at 100 GHz, as 
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Table 12. Total intensity and polarized flux densities, polarization fraction, and polarization position angles for Tau A measured 
from the Planck full-mission maps at 30 to 353 GHz, WMAP at 23 to 94 GHz (Weiland et al. 2011), and with the IRAM 30 m 
telescope at 89 GHz (Aumont et al. 2010). 


Freq 

I 

p 

PH 

Pol. Angle" 

Pol. Angle*’ 

Pol. Angle"’'* 

Pol. Angle^ 

[GHz] 

[Jy] 

[Jy] 

[%] 

[deg] 

[deg] 

[deg] 



30 

344.23 

+ 

0.27 

24.44 

+ 

1.05 

7.10 

+ 

0.33 

-84.54 

+ 

0.54 

+ 

0.50 

-83.71 

+ 

1.40 

+ 

0.50 

-89.26 

+ 

0.25 

+ 

0.50 






44 

292.68 

+ 

0.23 

19.07 

+ 

1.10 

6.51 

+ 

0.51 

-88.34 

+ 

0.32 

+ 

0.50 

-86.93 

+ 

0.47 

+ 

0.50 

-88.65 

+ 

0.79 

+ 

0.50 






70 

259.99 

+ 

0.11 

20.55 

+ 

0.61 

7.90 

+ 

0.32 

-84.24 

+ 

0.23 

+ 

0.50 

-85.03 

+ 

1.32 

+ 

0.50 

-87.49 

+ 

1.33 

+ 

0.50 






100 

215.16 

+ 

0.06 

15.54 

+ 

0.14 

7.22 

+ 

0.06 

-88.53 

+ 

0.11 

+ 

0.62 

-87.52 

+ 

0.13 

+ 

0.62 

-87.59 

+ 

0.26 

+ 

0.62 

-87.74 

+ 

0.26 

+ 

0.62 

143 

167.10 

+ 

0.04 

12.02 

+ 

0.08 

7.19 

+ 

0.05 

-84.85 

+ 

0.13 

+ 

0.62 

-85.72 

+ 

0.15 

+ 

0.62 

-87.03 

+ 

0.35 

+ 

0.62 

-87.22 

+ 

0.34 

+ 

0.62 

217 

124.21 

+ 

0.04 

10.09 

+ 

0.08 

8.12 

+ 

0.06 

-87.33 

+ 

0.12 

+ 

0.62 

-88.73 

+ 

0.18 

+ 

0.62 

-88.84 

+ 

0.55 

+ 

0.62 

-88.74 

+ 

0.55 

+ 

0.62 

353 

82.17 

+ 

0.67 

9.88 

+ 

0.17 

12.02 

+ 

0.23 

-86.11 

+ 

0.37 

+ 

0.62 

-85.15 

+ 

0.37 

+ 

0.62 

-85.16 

+ 

1.93 

+ 

0.62 

-88.38 

+ 

2.06 

+ 

0.62 

23 

383.80 

+ 

9.60 

27.17 

+ 

0.68 

7.08 

+ 

0.25 

-88.50 

+ 

0.10 

+ 

1.50 
















33 

342.80 

+ 

6.40 

23.80 

+ 

0.44 

6.94 

+ 

0.18 

-87.70 

+ 

0.10 

+ 

1.50 
















41 

317.70 

+ 

8.60 

22.12 

+ 

0.60 

6.97 

+ 

0.27 

-87.30 

+ 

0.20 

+ 

1.50 
















61 

276.00 

+ 

5.20 

19.31 

+ 

0.36 

7.00 

+ 

0.19 

-87.70 

+ 

0.40 

+ 

1.50 
















94 

232.80 

+ 

9.70 

16.60 

+ 

0.73 

7.13 

+ 

0.43 

-88.70 

+ 

0.70 

+ 

1.50 
















89 

195.00 

+ 

11.0 

14.50 

+ 

3.20 

8.80 

+ 

0.02 

-88.20 

+ 

0.20 

+ 

0.50 

















^ Position angle for Planck calculated using the maximum likelihood filtering method. 

^ Position angle for Planck calculated using aperture photometry. 

^ Position angle for Planck LFI channels (30, 44, and 70 GHz) calculated using the maximum likelihood filtering method on the special beam 
deconvolved maps. 

^ Position angle for Planck HFI channels (100, 143, 217, and 353 GHz) calculated using wider aperture photometry to allow for the angular extent 
of Tau A. 

^ Position angle for Planck HFI channels (100,143, 217, and 353 GHz) calculated using wider aperture photometry, and using HFI maps for which 
the diffuse bandpass correction has not been applied. 

Note: In the Planck channels, the statistical error bars in the polarization position angle of Tau A do not reflect the true uncertainties of these 
measurements. In addition to the statistical error, a 0^5 systematic error has to be added to the LFI measurements (Planck Collaboration III 2016) 
and 0°62 to the HFI measurements (Planck Collaboration VIII 2016) as indicated by the second + symbol. Similarly, the WMAP errors in the 
polarized position angle are statistical and a systematic error of 1°5 has been added (Weiland et al. 2011). 


compared with 561 + 156 mJy and 418 + 118mJy measured by 
PdBI, for PKS 0851-f202 and 3C273, respectively. 

We have also cross-matched the recent IRAM polarimet- 
ric survey from Agudo et al. (2014) at 86 GHz, again with the 
Planck 70 and 100 GHz PCCS2 catalogues. Although there are 
130 and 133 sources common to both samples, if we restrict the 
comparison to those with measured Planck polarized flux den¬ 
sities (rather than upper limits), we are left with five and seven 
sources at 70 and 100 GHz, respectively. At 70 GHz, two of the 
five sources (3C273 and 3C279) have similar polarized flux 
densities in both data sets, 559 and 522 mJy in IRAM, as com¬ 
pared with 519 + 93 and 368 ± 85 in Planck. At 100 GHz, three 
of the seven sources (3C273, 3C279, and PKS 1055+01) have 
similar polarized flux densities, 522, 559 and 305 mJy in IRAM 
as compared with 566+38,430+37, and 349 + 31 mJy in Planck, 
respectively. 

Additionally, we have compared the polarized flux densities 
found for 3C273, the only bright source in polarization in the 
sample of sources observed simultaneously with the VLA and 
Planck in the spring of 2013. The Planck 30 GHz polarized flux 
density (colour-corrected for comparison with the VLA) was 
854 + 82 mJy, while the VLA observed 843 + 50 mJy. At 44 GHz 
the Planck polarized flux density was 567 + 131 mJy, as com¬ 
pared with 623 + 70mJy seen by the VLA. 

3.6. Summary of validation 

The several internal and external validation tests described in 
Sect. 3 allow us to assess reliability and completeness as a func¬ 
tion of flux density in each Planck band. These tests also allow 
us to assess the accuracy of positions and flux densities tabulated 


in the PCCS2 and PCCS2E catalogues. With the possible excep¬ 
tion of a 6% or 1-1.5 cr difference between Planck flux densities 
measured in the noisy 44 GHz map and those measured from the 
ground, no clearly significant discrepancies in any of these quan¬ 
tities was found. We thus conclude that both positions and flux 
densities in the PCCS2 are valid within the tabulated statistical 
errors in total intensity. Note, however, that we have no direct test 
of the flux densities scales at 353 and 545 GHz. Regarding the 
measurements in polarization, the number of external surveys 
available at these frequencies in polarization is very limited. In 
particular, we have compared the flux densities in polarization of 
Tau A, the Crab nebula with recent measurements of WMAP and 
other high resolution instruments on the ground and we have not 
found significant discrepancies in the polarized flux densities. 
However, we And a small discrepancy in the measurement of the 
polarization angle of Tau A in some of the Planck channels. For 
this object, one of the brightest compact sources in polarization, 
there is a small amount of signal in the U maps at the position of 
the source, where little or no signal is expected for this object. 
This signal, much smaller than the signal in the Q map, does 
not have an effect in the measurement of the total polarized flux 
density since the Q and U flux densities are added in quadrature 
(as shown in Eq. 1), but it can explain the discrepancy in the 
polarization angle with respect to external measurements. Since 
our statistical errors do not account for this small systematic ef¬ 
fect, when we propagate the errors in the measurements of the 
polarized flux density into the errors in the measurement of the 
polarization angle, the errors that we obtain can be underesti¬ 
mated. 
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4. Characteristics of the PCCS2 

In Fig. 1, we displayed the sensitivity of the PCCS2 compared 
with the PCCS, the ERCSC, and several other CMB projects. For 
the PCCS2 we define the sensitivity to be the fiux density at the 
90 % completeness limit for each Planck channel. The improve¬ 
ments between the PCCS2 and the PCCS are most apparent for 
the LFI channels. This is to be expected given the larger increase 
of data for LFI than for HFI in the full mission. Additionally, for 
the higher frequency channels the foregrounds are a significant 
noise source for the detection of compact sources; the reduc¬ 
tion of the instrumental noise resulting from longer integration 
may therefore not increase the depth of the catalogue as much as 
might be expected. Finally, the estimated sensitivity of the cata¬ 
logue is worse at 857 GHz (Sect. 3.2.3), owing to the improved 
understanding of the reliability. 

Table 13 compares the characteristics of the PCCS2, the 
PCCS2E, and the PCCS. The total number of sources in each 
catalogue is given as well as the number outside the Galactic 
region. The numbers of sources in the extragalactic zone in gen¬ 
eral increase between the PCCS and the union of the PCCS2 and 
PCCS2E. For the highest three frequency channels, the union of 
the PCCS2 and PCCS2E catalogues contains many more sources 
than the PCCS due to the lower S/N threshold applied in the 
PCCS2E than in the Galactic region of the PCCS. In the extra- 
galactic zone, we compare the average uncertainty on the fiux 
density and the 90 % completeness values. We can see that the 
PCCS2 has lower uncertainties and is consequently more com¬ 
plete than the PCCS, except at 857 GHz, where the completeness 
has dropped (see Sect. 3.2.3). 

Table 14 shows the characteristics of the subset of sources 
with significant polarized emission, for both the PCCS2 and 
PCCS2E catalogues. The majority of the significantly polarized 
sources are in the Galactic plane region; hence for the HFI chan¬ 
nels the majority of these sources are in the PCCS2E catalogue. 

Table 15 shows the numbers of sources internally matched in 
adjacent frequency channels, within the union of the PCCS2 and 
PCCS2E. It shows the number of sources matched both above 
and below in frequency (e.g. sources at 100 GHz found in both 
the 70 and 143 GHz catalogues), those matched either above 
or below in frequency (a less stringent criterion), and the per¬ 
centage of sources so matched. Note that sources matched by 
the “above or below criterion” will include as a subset those 
sources meeting the more stringent “above and below” crite¬ 
rion. A source is considered to be matched if the positions are 
closer than the larger FWHM of the two channels. A catalogue 
was extracted from the IRIS 100 yum map (Miville-Deschenes & 
Lagache 2005) using the MHW2 pipeline, and used as the neigh¬ 
bouring channel above 857 GHz. The IRIS mask, which removes 
around 2.1 % of the sky, was applied to the 857 GHz catalogue 
before performing this comparison, and this reduces the number 
of sources in the union to 47 156, a decrease of about 2.1 %. 
The number of matches given for the 857 GHz channel only 
includes sources outside the IRIS mask. For the 30 GHz chan¬ 
nel, the matches were evaluated using only the channel above, 
44 GHz. The low percentage of internal matches of the 30 GHz 
channel results from two factors: the generally negative spectral 
index of the sources at these frequencies; and the relatively low 
sensitivity of the 44 GHz receivers. 

Figure 25 shows histograms of the spectral indices obtained 
via the non-blind fiux density extraction from the neighbour¬ 
ing channels. As expected, the high-frequency channels (545 
and 857 GHz) are dominated by dusty galaxies and the low- 
frequency ones are dominated by synchrotron sources, where the 



-10 12 3 


Spectral index 


Fig. 25. Histograms of spectral index for sources in PCCS2. The 
changes in the source populations with frequency are clearly vis¬ 
ible, the lower frequencies are dominated by synchrotron sources 
and the higher frequencies by dusty ones. At the intermediate 
frequencies both source populations are discernible. Between the 
top two panels, there is a visible shift in the peak of the his¬ 
togram. This is due to a steepening of the spectral indices of the 
radio sources. In each panel we give A, the number of sources in 
the histogram. 
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Fig. 26. Colour-colour plots. Red crosses represent sources from 
the PCCS2 and blue diamonds sources from the PCCS2E. Top: 
common frequency 217 GHz. We can see the non-thermal and 
thermal source populations of the PCCS2; the PCCS2E con¬ 
tains significantly more thermal sources than the PCCS2, as ex¬ 
pected given the that PCCS2E contains the Galactic plane re¬ 
gion. Bottom: common frequency 857 GHz. The PCCS2 is con¬ 
sistent with a population of cold sources spanning a narrow 
range in temperature, whereas the PCCS2E shows a wider dis¬ 
tribution of source properties. 


change in the dominant source population occurs between 217 
and 353 GHz. It can also be seen that there is a shift in the peak of 
the histogram between the top two panels, between 30-44 GHz 
and 44-70 GHz and above. The reason for this is a steepening of 
the spectral indices of radio sources, which has been seen pre¬ 


viously (Massardi et al. 2009; Planck Collaboration XIII 2011; 
Planck Collaboration XXVIII 2014). 

Eigure 26 shows colour-colour plots for sources from the 
PCCS2 and the PCCS2E. The positions in the catalogue at the 
common frequency in the colour-colour plot are used to perform 
a non-blind extraction of the flux densities from the maps at the 
neighbouring frequencies. These are used together with the com¬ 
mon frequency-channel flux densities to construct the colour- 
colour plot. The common frequency in the top panel is 217 GHz, 
and here we can see the two populations, which are also seen 
in the spectral-index histogram (Eig. 25), for the PCCS2. Also 
apparent is the domination of the PCCS2E at this frequency 
by thermal sources, since the PCCS2E includes sources in the 
Galactic plane. The bottom panel uses the IRIS map and the 
Planck 545 GHz maps for the non-blind extraction of flux den¬ 
sities at the locations of the 857 GHz catalogues. The PCCS2 
colour-colour distribution is consistent with a population of cold 
sources with a narrow range of temperatures (around 10-40 K) 
in the range of spectral indices is shown in Eig. 25. The larger 
dispersion of values in the PCCS2E by comparison, is indicative 
of the greater noise levels in the region of sky corresponding to 
the PCCS2E as well as a broader distribution of temperatures 
for the sources. 


5. The PCCS2: access, content and usage 

The PCCS2 is available from the Planck Legacy Archive.^ It 
is composed of 15 single frequency catalogue EITS files, one 
per LEI channel and two per HEI channel. In addition there are 
associated maps, again provided as EITS files, which are de¬ 
scribed further in Sect. 5.1. Additional information about the 
catalogue content and format can be found in the Explanatory 
Supplement,^ in the EITS file headers, and in the first PCCS pa¬ 
per (Planck Collaboration XXVIII 2014). Here we summarize 
the catalogue contents, focusing on the additional features of the 
PCCS2 catalogues. 

- Source identification: NAME (e.g. PCCS2 

030G184.54-05.78). 

- Position: GLON and GLAT contain the Galactic coordinates, 
and RA and DEC give the same information in equatorial 
coordinates (J2000). 

- Elux density: the four estimates of flux density (DETELUX, 
APERELUX, PSEELUX, and GAUELUX) in mJy, and their 
associated uncertainties. 

- Source shape: the elliptical Gaussian fit to the source; i.e. the 
semi-axes, and orientation. 

- Polarization measurements: the polarized fiux density and 
polarization angle and their associated errors for signifi¬ 
cantly polarized sources; provided for all seven of the nine 
Planck channels that have polarization data (30-353 GHz). 
The polarization angles are defined as increasing anticlock¬ 
wise (north through east) following the lAU convention; the 
position angle zero is the direction of the north Galactic pole. 

- Marginal polarization measurements: measurements for less 
significantly polarized sources, as described in Sect. C.2; 
these are provided for the 100-353 GHz channels only. 

- Source extension: the EXTENDED fiag is set to 1 if a source 
is extended. A source is extended if the geometric mean of 
the elliptic Gaussian fit to the source is greater than one-and- 
a-half times the fitted EWHM from Table 2. 

^ http://WWW.cosmos.esa.int/web/planck/pla 

^ http://wiki.cosmos.esa.int/planckpla2015 
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Table 13. PCCS2 and PCCS2E characteristics compared with those of the PCCS. 


Channel 

30 

44 

70 

100 

143 

217 

353 

545 

857 

Freq[GHz]. 

28.4 

44.1 

70.4 

100.0 

143.0 

217.0 

353.0 

545.0 

857.0 

A [fim] . 10561 

6807 

4260 

3000 

2098 

1382 

850 

550 

350 

Number of sources 










PCCS2 . 

1560 

934 

1296 

1742 

2160 

2135 

1344 

1694 

4891 

PCCS2E. 




2487 

4139 

16842 

22665 

31068 

43290 

Union PCCS2+PCCS2E. 




4229 

6299 

18977 

24009 

32762 

48181 

PCCS^ . 

1256 

731 

939 

3850 

5675 

16070 

13613 

16933 

24381 

Number of sources 










- in extragalactic zone^ 










PCCS2 . 

745 

367 

504 

1742 

2160 

2135 

1344 

1694 

4891 

PCCS2E. 




0 

0 

26 

289 

839 

2097 

Union PCCS2+PCCS2E. 




1742 

2160 

2161 

1633 

2533 

6988 

PCCS^ . 

572 

258 

332 

1483 

1779 

1745 

1424 

3566 

7270 

Flux densities [mJy] 










- in extragalactic zone^ 










PCCS2 : minimum^. 

376 

603 

444 

232 

147 

127 

242 

535 

720 

: 90% completeness . . 

426 

676 

489 

269 

177 

152 

304 

555 

791 

: uncertainty. 

87 

134 

101 

55 

35 

29 

55 

105 

168 

PCCS2E: minimum^. 






189 

350 

597 

939 

: 90% completeness . . 






144 

311 

557 

927 

: uncertainty. 






35 

73 

144 

278 

PCCS^ : minimum^. 

461 

825 

566 

267 

169 

140 

273 

445 

668 

: 90% completeness . . 

575 

1047 

776 

300 

190 

180 

330 

570 

680 

: uncertainty. 

109 

198 

149 

62 

39 

33 

65 

119 

188 

^ Planck Collaboration XXVIII (2014) 










^ 30-70 GHz: as in PCCS, the extragalactic zone is 

given by \b\ > 30°. 100-857 GHz: outside of Galactic region where the reliability cannot be 

accurately assessed. Note that for the PCCS2E the only sources that occur in this region lie in the filament mask. 



^ Minimum fiux density of the catalogue 

in the extragalactic 

zone after excluding the faintest 10 % of 

' sources. 




Table 14. PCCS2 and PCCS2E polarization characteristics for sources with polarized emission with significance > 99.99%. 


Channel 



30 

44 

70 

100 

143 

217 

353 

Number of significantly polarized sources in PCCS2 


122 

30 

34 

20 

25 

11 

1 

Minimum polarized flux density^ [mJy] . 



117 

181 

284 

138 

148 

166 

453 

Polarized flux density uncertainty [mJy]. 



46 

88 

91 

30 

26 

30 

81 

Minimum polarized flux density completeness 90% [mJy] . . 

199 

412 

397 

135 

100 

136 

347 

Minimum polarized flux density completeness 95% [mJy] . . 

251 

468 

454 

160 

111 

153 

399 

Minimum polarized flux density completeness 100% [mJy] . 

600 

700 

700 

250 

147 

257 

426 

Number of significantly polarized sources in PCCS2E. 




43 

111 

325 

666 

Minimum polarized flux density^ [mJy] . 






121 

87 

114 

348 

Polarized flux density uncertainty [mJy] . 






52 

44 

55 

178 

Minimum polarized flux density completeness 90% [mJy] . . 




410 

613 

270 

567 

Minimum polarized flux density completeness 95% [mJy] . . 




599 

893 

464 

590 

Minimum polarized flux density completeness 100% [mJy] . 




835 

893 

786 

958 


^ Minimum polarized flux density of the catalogue of significantly polarized sources after excluding the faintest 10 % of sources. For the LFI 
channels we have not considered the sources that have been flagged as unidentified. There are nine, one, and one of these unidentified sources at 
30, 44, and 70 GHz, respectively. 


- External validation: EXT_VAL contains a summary of the 
inter-channel and external validation. See the definition be¬ 
low. 

“ Positional coincidence identification with a previous Planck 
catalogue: the ERCSC and PCCS columns indicate the 


names of the ERCSC and PCCS counterparts, if they exist, 
at that channel. 

- Degree of reliability: in the PCCS2 catalogue the 
HIGHEST_RELIABILITY_CAT column contains the high¬ 
est reliability catalogue to which the source belongs. 
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Table 15. Sources matched between neighbouring channels. 
Given that, for the HFI channels, the sky area corresponding to 
the PCCS2 and PCCS2E is different for every channel, the com¬ 
parison between frequency channels must necessarily be per¬ 
formed on the union of the PCCS2 and PCCS2E. The fraction 
matched from the same analysis applied to the previous PCCS is 
shown in brackets in the last column, for comparison purposes. 
Note that 217 GHz is different from other bands in the ratio of 
the number of above-and-below matches compared to the num¬ 
ber of above-or-below. This is of course because of the change of 
spectral index at that point, as discussed in Sect. 4 and Figs. 25 
and 26. 


No. matched Fraction matched 


No. 

Channel sources 

Above 

and 

below 

Above 

or 

below 

PCCS2 

and 

PCCS2E 

(PCCS) 

30V . . 

1560 


799 

51.2% 

(50.1%) 

44 ... . 

934 

700 

851 

91.1% 

(90.8%) 

70 ... . 

1296 

735 

1113 

85.9% 

(86.8%) 

100 ... . 

4229 

1047 

3049 

72.1% 

(71.6%) 

143 ... . 

6299 

2734 

5163 

81.9% 

(81.9%) 

217 ... . 

18977 

3837 

14928 

78.7% 

(66.1%) 

353 ... . 

24009 

12171 

20867 

86.9% 

(88.7%) 

545 ... . 

32762 

17003 

28423 

86.8% 

(85.8%) 

857^ . .. 

47156 

14578 

35390 

75.0% 

(74.9%) 


^ The 30 GHz channel is only matched with the 44 GHz channel above. 
^ The 857 GHz channel is matched above with a catalogue extracted 
from the IRIS maps using the HFI-MHW. Both catalogues are cut 
with the IRIS mask prior to matching. 


“ Reason for inclusion in PCCS2E: the WHICH_ZONE flag 
encodes why the source has been placed in the PCCS2E. 

- Cirrus indicators: a fraction of the sources detected in the up¬ 
per HFI bands could be associated with Galactic interstellar 
medium features or cirrus. For the 217-857 GHz channels, 
the CIRRUS_N and the new column, SKY_BRIGHTNESS, 
may be used as cirrus indicators. The CIRRUS_N column 
is defined as in PCCS, the number of sources detected at 
857 GHz (using a uniform S/N threshold of 5) within a 1° ra¬ 
dius of the source. The SKY .BRIGHTNESS is defined as the 
mean 857 GHz brightness within a 2° radius of the source. 
See Sect. 5.2 for details. 

Note that all flags are evaluated per frequency channel. Two 
flags require information from the 857 GHz channel for their 
evaluation, CIRRUSJ^ and SKY.BRIGHTNESS. These are, 
however, evaluated independently for each frequency channel 
for which they are provided. This means, for example, a source 
that as been observed at 545 and 857 GHz at slightly different po¬ 
sitions in each channel (but close enough to be considered to be 
the same source), may have different values for the CIRRUS_N 
and SKY_BRIGHTNESS flags in each channel. 

The EXT.VAL column summarizes the cross-identification 
with external catalogues. Its definition has been modified 
slightly with respect to the PCCS. The EXT.VAL flag now has 
a value of 0, 1, 2, or 3 as described below: 

0: The source has no clear counterpart in any of the external 
catalogues and it has not been detected in a neighbouring 
Planck channel. 


1: The source has no clear counterpart in any of the external 
catalogues, has not been detected in a neighbouring Planck 
channel, but was detected in the same channel in the PCCS. 
2: The source has no clear counterpart in any of the external 
catalogues, but it has been detected in a neighbouring Planck 
channel in this release. For the HFI channels, we consider the 
catalogues extracted from the IRIS maps as a neighbouring 
Planck channel, given the common detection algorithm ap¬ 
plied to both data sets. 

3: The source has a clear counterpart in one of the radio cat¬ 
alogues (CRATES, Healey et al. 2007; NEWPS, Lopez- 
Caniego et al. 2007; AT20G, Murphy et al. 2010), the 
Revised IRAS-FSC Redshift Catalogue (RIFSCz; Wang 
et al. 2014), or the submillimetre catalogue of H-ATLAS 
(Bales et al. 2010). 

This flag provides extra information about the reliability of in¬ 
dividual sources: those flagged as EXT.VAL = 3 are already 
known; those with EXT.VAL = 2 (or 1) have been detected in 
other Planck channels (or maps) and are therefore potentially 
new sources, and those with EXT.VAL = 0 appear in only a sin¬ 
gle channel and only in this release, and hence are more likely 
to be spurious. 

The EXT.VAL flag requires data from the neighbouring 
channels for its evaluation. It is, however, evaluated per fre¬ 
quency channel. If, for instance, a source is identified with an 
external catalogue and a source in a neighbouring channel, the 
source in the neighbouring channel will only be identified with 
the external catalogue if it also satisfies the identification criteria 
with the external catalogue. 

It should be noted that there is no column that contains the 
coordinate uncertainties for each source. The errors in position 
are purely statistical and may be determined for each source us¬ 
ing Eq. (7) with the parameters given in Table 8. There is also no 
column that contains the S/N of the detection, since this is given 
by DETFLUX/DETFLUX_ERR. 

5.1. Maps associated with the catalogues 

Along with the source catalogues we provide associated maps 
for the HFI channels. There are three types of maps provided, 
which are shown in Fig 8. These are a zone mask, a noise level 
map, and a S/N threshold map. 

The zone map shows the areas of the sky covered by the 
PCCS2 and PCCS2E catalogues. The zone map takes the value 
of zero in the areas corresponding to the PCCS2 and is non¬ 
zero in areas corresponding to the PCCS2E. It is related to the 
WHICH_ZONE flag, in that a non-zero value encodes the reason 
why sources in that patch of sky are placed in the PCCS2E. The 
value 1 corresponds to the filament mask outside of the Galactic 
region, 2 corresponds to the Galactic region, and 3 corresponds 
to the filament mask inside the Galactic region. 

The noise level map corresponds to the detection noise for 
compact sources. This is not the same as the instrumental noise, 
because it includes “noise” from all signals other than compact 
sources. 

The S/N threshold map contains the thresholds required at 
each location on the sky to produce a catalogue of the stated 
reliability. Threshold maps for the 80 %, 85 %, 90 %, and 95 
% reliability catalogues are provided. The S/N threshold for the 
80% reliability map contains the S/N cut applied to the PCCS2E 
in the area of the sky where the zone map is non-zero. The S/N 
threshold maps for higher reliabilities contain null values in this 
region. 
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In total there are, for each HFI channel, six associated maps: 
one zone map, one noise level map, and four S/N threshold maps. 
As described in Sect. 3.1.2, the completeness of the PCCS2 and 
PCCS2E may be evaluated using these maps. By providing the 
S/N threshold maps for the higher reliabilities the complete¬ 
nesses for the higher reliability subsets of the PCCS2 may also 
be evaluated. Indeed the completeness for a given reliability and 
region of the sky may be assessed using these data. 

5.2. Cautionary notes on the use of catalogues 

The PCCS2 supersedes previous Planck compact source cata¬ 
logues (ERCSC and PCCS), because it has been produced us¬ 
ing the full mission data and the latest processing and calibra¬ 
tion pipelines. Since the three sets of catalogues have been pro¬ 
duced from the analysis of maps that average different amounts 
of data, some idea of the variability of the sources could be ob¬ 
tained from the comparison of the three, although for this pur¬ 
pose it would be better to analyse the single survey maps, the 
time-ordered data or the user specified time interval maps that 
can be obtained from the Planck Legacy Archive. 

As noted earlier, the aim of the PCCS2 is to provide as com¬ 
plete a list as possible of Planck sources with a reasonable, and 
user-adjustable degree of reliability. The criteria used to include 
or exclude candidate sources differ from channel to channel and 
in different parts of the sky; they also are based on different S/N 
levels between channels and as a function of position on the sky. 
These differences are consequences of our desire to make the 
catalogue as complete as possible, yet maintain > 80 % reliabil¬ 
ity; the differences have to be taken into account when using the 
PCCS2 for statistical studies. 

We now turn to several specific cautions and comments for 
users of the PCCS2 and PCCS2E. 

Bandpass corrections to polarization : Eor many sources in the 
three lowest Planck frequency channels, the bandpass correction 
of the Q and U flux densities is not negligible. Even though we 
have attempted to correct for this effect on a source by source 
basis and have propagated this uncertainty into the error bars on 
the polarized flux densities and polarization angles, there is still 
room for improvement. This can be seen in the residual leakage 
present at the position of Tau A in the Stokes U maps. It is an¬ 
ticipated that there will be future updates to the LEI catalogues 
once the bandpass corrections and errors have been improved. 

Variability: At radio frequencies, up to and including 217 GHz, 
many of the extragalactic sources are variable. The measure¬ 
ments of their flux densities provided in the catalogues are, how¬ 
ever, averages over the full Planck mission. It should be noted 
therefore that follow-up observations of these sources may show 
significant differences from those provided. 

Contamination from CO: At infrared/submillimetre frequen¬ 
cies (100 GHz and above), the Planck bandpasses straddle ener¬ 
getically significant CO lines. The effect is the most significant 
at 100 GHz, where the line might contribute more than 50 % of 
the measured flux density for some Galactic sources. 

Photometry: Each source has multiple estimates of flux den¬ 
sity: DETFLUX, APERFLUX, GAUFLUX, and PSFFLUX, as 
defined above. The evaluation of APERFLUX makes the small¬ 


est number of assumptions about the data and hence is the most 
robust, especially in regions of high non-Gaussian background 
emission, but it may have larger uncertainties than the other 
methods. Hence, a general recommendation for which estimate 
to use for unresolved sources would be DETFLUX for 30 to 
217 GHz and APERFLUX for 353 to 857 GHz. Note that for 
a specific source the nature of the local background will influ¬ 
ence the best choice for the flux estimator. For bright resolved 
sources, GAUFLUX is recommended, taking int account that it 
may not be robust for sources close to the Galactic plane due 
to the strong backgrounds. In the PCCS2 and PCCS2E we pro¬ 
vide polarized flux densities for two methods, one obtained from 
the measured flux densities in the filtered maps of Q and U and 
the other obtained from the measured flux densities on the unfil¬ 
tered maps of Q and U using aperture photometry. Both meth¬ 
ods agree for the brightest sources but, because the noise level 
is higher in the unfiltered maps, for weaker sources the polar¬ 
ized flux densities obtained with the filtering method are more 
robust. In addition, we have found that at the position of very 
bright sources in polarization one can see a spurious signal in¬ 
troduced by the complex beams in polarization. This spurious 
signal is small compared with the flux density of the sources, 
but in cases like Tau A, where most of the signal is in the Q map, 
this signal can have an impact on the flux density measured in the 
U map, which is particularly important when calculating the po¬ 
larization position angle. For this reason position angles should 
be used with caution. 

Calibration: The absolute calibration uncertainties of Planck 
are well below 1 % for 30-353 GHz (Planck Collaboration V 
2016; Planck Collaboration VIII 2016), while for 545 and 
857 GHz the absolute calibration uncertainty is < 7 %, which 
is primarily due to a 5 % systematic uncertainty arising from the 
planet models (Planck Collaboration VIII2016). This systematic 
uncertainty is not included in the internal validation (not simu¬ 
lated) or in the comparison with Herschel data (which use the 
same planet models; Griffin et al. 2010). In addition, there is a 
0?5 and 0?3 systematic uncertainty in the polarization position 
angles in the LFI and HFI, respectively. 

Colour-correction: The flux-density estimates have not been 
colour-corrected because this implies fixing the spectral index 
that describes each source, and the user may want to apply 
a colour-correction based on a specific spectral index deter¬ 
mined with a higher-resolution experiment. In the Planck Legacy 
Archive there is a tool to apply the colour-correction to the 
source flux density for a user supplied spectral index. Colour- 
corrections are described in Planck Collaboration II (2016) and 
Planck Collaboration VII (2016). Note that the term bandpass 
correction in this paper refers to the correction required due to 
the mismatch in the bandpass between orthogonally polarized 
detectors and not to the colour-correction of the flux density in 
Stokes /. 

Cirrus/ISM: The upper bands of HFI could be contami¬ 
nated by apparent sources associated with Galactic interstellar 
medium features (ISM) or cirrus. The values of the parameters 
CIRRUS_N and SKY_BRIGHTNESS can be used as indicators 
of contamination. CIRRUS JN can be used to flag sources that 
might be clustered together and thereby associated with ISM 
structure. In order to provide some indications of the range of 
values of these parameters that could indicate contamination, we 
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compared the properties of the IRAS-identified and non-IRAS- 
identified sources for both the PCCS2 and the PCCS2E, outside 
the Galactic plane. At Galactic latitudes \b\ > 20°, we can use 
the RIFSCz (Wang et al. 2014) to provide a guide to the likely 
nature of sources. 

We compare the PCCS2 857 GHz catalogue and the PCCS2E 
857 GHz catalogue with the IRAS sources in the RIESCz us¬ 
ing a Sarcmin matching radius. Of the 4891 sources in the 
PCCS2 857 GHz catalogue, 3094 have plausible IRAS coun¬ 
terparts, while 1797 do not. Examination of histograms of the 
CIRRUS_N and SKY_BRIGHTNESS parameters in the PCCS2 
show that these two classes of objects behave rather differ¬ 
ently. The IRAS-identified sources have a peak sky brightness 
at about 1 MJy sr“^. The non-IRAS-identified sources have a bi- 
modal distribution with a slight peak at 1 MJy sr“^ and a sec¬ 
ond peak at about 2.6MJysr“^ Both distributions have a long 
tail, but the non-IR AS-identified tail is much longer. On this ba¬ 
sis, sources with SKY_BRIGHTNESS > 4MJysr“^ should be 
treated with caution. In contrast non-IR AS-identified sources 
with SKY_BRIGHTNESS < 1.4 MJy sr”! are likely to be reli¬ 
able. Examination of their sky distribution, for example, shows 
that many such sources lie in the IRAS coverage gaps. The 
CIRRUS_N flag tells a similar story. Both IRAS-identified and 
IRAS non-identified sources have a peak CIRRUS _N value of 
2, but the non-identified sources have a far longer tail. Very few 
IRAS-identified sources have a value > 8 but many unidentified 
sources do. These should be treated with caution. 

The PCCS2E 857 GHz catalogue contains many more 
sources with \b\ > 20° of which 1235 are identified with IRAS 
sources in the RIFSCz and 9235 are not. As with the PCCS2 cat¬ 
alogue the distributions of CIRRUS _N and SKY_BRIGHTNESS 
are different, where the differences are even more pronounced 
for these PCCS2E sources. Once again, few IRAS-identified 
sources have SKY_BRIGHTNESS > 4 MJy sr“\ but the uniden¬ 
tified sources have brightnesses extending to >55MJysr“^ 
Similarly, hardly any of the IRAS-identified sources have 
CIRRUS_N > 8, but nearly half the unmatched sources do. Of 
the 9235 PCCS2E 857 GHz sources that are not identified with 
an IRAS source and that lie in the region \b\ > 20°, 1850 (20 %) 
have WHICH_ZONE = 1, 2637 (29 %) have WHICH_ZONE = 
2, and 4748 (51 %) have WHICH_ZONE = 3. The PCCS2E cov¬ 
ers 30.36 % of the region \b\ > 20°, where 2.47 % is in the fila¬ 
ment mask, 23.15 % in the Galactic region and 4.74 % in both. 
If the 9235 unidentified detections were distributed uniformly 
over the region \b\ > 20°, we can predict the number of uniden¬ 
tified sources in each zone and compare this to the values we 
have. We find that there are 2.5 and 3.3 times more sources than 
expected in zones 1 and 3, showing that the filament mask is in¬ 
deed a useful criterion for regarding sources detected within it 
as suspicious. 

It should be noted that the EXTENDED flag could also be 
used to identify ISM features, but nearby Galactic and extra- 
galactic sources that are extended at Planck spatial resolution 
will also meet this criterion. 

6. Conclusions 

The Second Planck Catalogue of Compact Sources has been pro¬ 
duced using the Planck full mission data. The catalogue lists 
sources detected in total intensity in each of its nine frequency 
bands between 30 and 857 GHz and polarization measurements 
at the positions of these sources for seven of the frequencies 
between 30 and 353 GHz. Its format has changed with respect 
to the ERCSC and the PCCS. We have divided the catalogue 


into two parts, the PCCS2 and PCCS2E, based on our ability 
to provide a measure of the reliability of each source detected 
at 100 GHz and above, where the available external catalogues 
of compact sources are not able to fully assess the reliability of 
the detections. Sources located inside the defined Galactic plane 
masks or situated along dusty filamentary structures (as defined 
in the cirrus masks) are in the PCCS2Ebecause the uncertainties 
in the number counts of the Galactic sources and the difficulty of 
simulating the diffuse dust emission near the beam scale in the 
higher frequency channels do not allow us to achieve the nec¬ 
essary consistency between the catalogues of input and detected 
sources in our reliability assessments. 

Given the increase in the volume of data between the nom¬ 
inal mission and the full mission, and the improvements in the 
data processing and calibration of the frequency channel maps, 
the PCCS2 supersedes previous Planck catalogues. The new cat¬ 
alogue is more complete than the PCCS, in particular for the LEI 
channels, owing to the large increase in the data available, eight 
sky surveys compared with the two and a half sky surveys of 
PCCS. In addition, improvements have been made in some of 
the techniques used to perform the photometry analysis, and in 
the reliability assessment of the catalogues. The completeness 
of the 857 GHz channel, however, has not improved because 
improvements in the reliability assessment resulted in a higher 
S/N threshold being applied in the formation of this catalogue. 
It should be noted, however, that the quality of the PCCS2 cata¬ 
logue at this channel is better than that of the earlier PCCS due 
to its greater reliability. 

The division of the HFI catalogues into the PCCS2 and 
PCCS2E has permitted the addition of a parameter in the PCCS2 
catalogue that will allow a user to define subsets of the cata¬ 
logue with higher reliability levels than the target integral reli¬ 
ability of 80%. Associated maps are provided that will allow 
the user to evaluate the completeness of their chosen reliabil¬ 
ity subset, or indeed of the catalogue as a whole. This added 
functionality gives the users of the PCCS2 the option of extract¬ 
ing high-reliability subcatalogues, and, in addition, provides a 
much more complete full catalogue, allowing studies of more 
sources and to fainter fiux densities. There are ongoing efforts in 
the Planck Collaboration to produce multifrequency catalogues 
that will complement the PCCS2; when completed, they will be 
available in the Planck Legacy Archive. 
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Appendix A: PSF photometry 

The flux density is obtained by fitting a model of the PSF at 
the position of the source. The model has four free parameters: 
the amplitude of the source, a background offset, and two co¬ 
ordinates for the location of the source. The PSF is obtained 
from the effective beam (Planck Collaboration II 2014; Planck 
Collaboration VI 2014). The model of the source is 

m = AP-\-C, (A.l) 

where P is the PSF at the position of the source (the integrated 
response to a point-like source), A is the amplitude of the source, 
and C is the (constant) background. The PSF at the position of 
the source is obtained from the effective beam, which is defined 
only at the centre of each map pixel, by means of a bicubic in¬ 
terpolation between adjacent pixels. This step is new: in the pre¬ 
vious version of the PCCS, the PSF was built from the effective 
beam at the centre of the pixel associated with the location of the 
source. The PSF model P depends therefore on the position of 
the source, P = P(Xs,ys)- 

The best-fit values of the parameters {3 = (A,C, Y^,y^) are 
found by minimizing the;^^ between the model and the data, d, 

VOS) = ~ {d - m), (A.2) 
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Fig.B.l. Comparison of injected and recovered MHW2 flux 
densities for point sources at 100 GHz. The saturation of points 
is proportional to the S/N of the MHW2 detection. The outliers 
towards the upper left are associated with faint sources with in¬ 
jected flux densities < 1 Jy located in areas with complicated 
backgrounds close to the Galactic plane. It is known that in these 
regions the MHW2 algorithm may give biased flux-density esti¬ 
mates. 

where N is the covariance matrix of the noise. The noise is as¬ 
sumed to be uncorrelated between pixels. The overall normal¬ 
ization of the noise is adjusted by setting = 1 at the best-flt 
value of the parameters. We also include the uncertainty of the 
background flux in the error of the flux density estimation. This 
has the effect of inflating the uncertainties to account for any 
mismatch between the modelled PSF and the true shape of the 
source and the background in the map. The uncertainties on the 
parameters are computed from the curvature of the;^^. The best- 
flt amplitude and its uncertainty are converted to units of flux 
density using the area of the PSF and the unit conversion from 
Kcmb to MJy sr“^ for each Planck channel. 

Appendix B: Gaussian fitting method 

To recover the shape of the source and its orientation, and to 
improve on the MHW2 flux-density estimate, Gaussian fitting 
is used. Taking point source parameters from the MHW2 cata¬ 
logue as input, a 2D Gaussian fit is performed at the location of 
each catalogue source. Six parameters are fitted: Galactic coor¬ 
dinates of the source (/, b)\ the flux density; the major and minor 
semi-axes; and an orientation angle. The source semi-axes and 
orientation angle are used to flag elongated sources. 

B.1. Downhill sinnplex in multidimensions 

The downhill simplex method in multidimensions, i.e. the 
Nelder-Mead or “amoeba” method (Press et al. 1992), which 
is useful for problems where the derivatives are not known, is 
used for optimization in multiparameter space. The functional to 
optimize is based on the reduced log-likelihood with some prior 
regularization for the size of the source defined by the effec¬ 
tive beam at the given frequency. The algorithm starts optimiza¬ 
tion at the MHW2 source location, assuming initially a circular 


Fig.B.2. Comparison of injected and recovered Gaussian flux 
densities for point sources at 100 GHz. The saturation of points 
is proportional to the S/N of the MHW2 detection. Gaussian 
fitting improves the outlier flux densities slightly, but the flux- 
density estimates are still biased. 
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Fig.B.3. Comparison of injected and recovered MHW2 flux 
densities for cold cores at 100 GHz. The saturation of points is 
proportional to the S/N of the MHW2 detection. The visible dig¬ 
itization is because only integer values were used for the flux 
densities of injected sources. 


Gaussian source with the FHWM of the PSF. Optimization of all 
six parameters converges in a reasonable number of iterations, 
usually less than 1000. 

B.2. Error estimation 

The downhill simplex method does not provide any information 
on the uncertainties in the optimized parameters. We estimate 
errors using a Markov Chain Monte Carlo (MCMC) method to 
sample the posterior distribution. In our case, the bum-in time 
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Fig.B.4. Comparison of injected and recovered Gaussian flux 
densities for cold cores at 100 GHz. The saturation of points is 
proportional to the S/N of the MHW2 detection. The visible dig¬ 
itization is because only integer values were used for the flux 
densities of injected sources. 


methods give almost the same point source flux-density values 
(Figs. B.l and B.2), but MHW2 gives less scatter for the fainter 
point sources. For the cold cores, Gaussian fitting works bet¬ 
ter and gives much less bias, especially for the bright sources 
(Figs. B.3 and B.4). The visible digitization of the flux densities 
reflects the fact that only integer values were used for the flux 
densities of injected sources. 

The group of outliers in the upper left part of 
Figs. B.l and B.2 correspond to faint point sources (< 1 Jy) 
in regions with complicated backgrounds, near the Galactic 
plane. Both MHW2 and Gaussian fitting fail to recover their 
flux densities accurately and give biased estimates. 

It is clear that Gaussian fitting is preferable for extended ob¬ 
jects like cold cores, giving more accurate flux densities as well 
as recovering the actual source shape and orientation. 


Appendix C: Bandpass mismatch and poiarization 
measures 

For sources with spectra differing from the CMB spectrum, 
bandpass mismatch causes leakage of temperature to polariza¬ 
tion (see Sect. 2.5.1). To correct for this leakage, we require a 
model of the source spectrum. That requirement was handled 
differently by LFI and HFI. 


can be very small since to start the chain the simplex-optimized 
value of the parameter from the previous step can be used, and 
it is usually very close to the peak of the distribution. The length 
of the MCMC chain is set to 1000 samples by default. 

To obtain the proper sampling in the MCMC chain it is nec¬ 
essary to define the optimal step size in the parameter space. 
This determination is complex since there are more than 100 000 
sources in the catalogue to process; it would be impossible to ad¬ 
just every MCMC chain, since it would take an excessive amount 
of time to try several values of the sampling step using an iter¬ 
ative approach. Hence we decided to use an analytical approxi¬ 
mation for the errors of most of the Gaussian profile parameters, 
following Hagen & Dereniak (2008). We used an analytical ap¬ 
proximation for the separable 2D Gaussian profile model; the er¬ 
rors are derived for five parameters; namely position coordinates, 
flux density, and two non-equal width values. The calculated er¬ 
rors were used to And the value of the step size in the parameter 
space to start the MCMC sampler, which delivered the error on 
the tilt angle. This combined technique allowed us to avoid most 
of the problems related to undersampled or oversampled MCMC 
chains. 

B.3. Validation of the method using model data 

The results of the Gaussian fit were compared with the MHW2 
positions and flux densities for a test set of sources with known 
parameters. The PR2 100 GHz map with injected sources and 
cold cores was used to check the 2D Gaussian fitting algorithm. 
The cold cores were modelled as 2D Gaussian profiles with uni¬ 
formly distributed random orientation angles. The major FWHM 
was in the range from 4'5 to 19', and the ellipticity parameter, 
e = FWHMmajor/FWHMminor. Varied from 1 (circular source) 
to 7 (highly elliptical source). All sources were uniformly dis¬ 
tributed across the sky. 

Plots comparing the flux density values for MHW2 and 
Gaussian fitting for this test set are shown in Figs. B.1-B.4. Both 


C.1. LFi-specific details 

In the case of LFI maps corrected for bandpass-mismatch leak¬ 
age, the spectral model is based on the diffuse component sep¬ 
aration analysis (Planck Collaboration II 2016). These mod¬ 
els are not particularly accurate for compact sources; more¬ 
over, the analysis is done at a common resolution of 1° FWHM, 
whereas for compact sources it is best to work at full resolution. 
Therefore the PCCS2 polarization measurements were evaluated 
by extracting the Q and U flux densities from the full-resolution 
uncorrected maps, and correcting for leakage at the flux density 
level, using 

= l§) -[pA{a-acM^)I, (Cl) 

corrected \ / raw \ ^ / 

where I is the total intensity, a = J In //J In v is the source spec¬ 
tral index in the relevant frequency band, and acuB is the spec¬ 
tral index of the CMB fluctuations (1.96,1.90, and 1.75, at 28.4, 
44.1, and 70.4 GHz, respectively). Pq^u is the projection factor 
derived in Planck Collaboration II (2016), which is evaluated at 
each map pixel. For the point-source correction we averaged the 
projection factors over the source pixels using the same weights 
as were used to extract the Q and U fluxes. 

C.1.1. Evaluation of LFI polarization measurements 

The polarized flux density P is calculated using Eq. (1), where 
the Q and U maps are obtained by applying the “Altered fusion” 
maximum-likelihood estimator (Argiieso et al. 2009; Lopez- 
Caniego et al. 2009) to the original Stokes Q and U maps. The 
errors in P are calculated using Eq. (4), adding in quadrature the 
error estimate obtained from the Monte Carlo simulations and 
calculating the 1 cr asymmetric errors for a Rayleigh distribution 
in the intervals [0,0.159] and [0.841,1]. The errors in the polar¬ 
ization angle, if/, are calculated using Eq. (5). 
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Fig. C.l. Comparison of polarized flux densities from FFP8 sim¬ 
ulations with (S'bpm) and without (^nobpm) bandpass mismatch 
leakage. The difference between the polarized flux densities 
normalized by the uncertainty on the polarized flux density, 
A = (^bpm - Snobpm)/o', is plotted against the input polarized 
flux density. The errors in the polarized flux density due to leak¬ 
age are always subdominant with respect to the uncertainties due 
to the noise. 


C.2. HFI-specific details 

Data from a number of bolometers are combined to make the 
HFI polarized frequency maps. Mismatch between the band- 
passes of the bolometers causes leakage from intensity to polar¬ 
ization for any source of emission which has a non-CMB spec¬ 
trum. Small uncertainties in the measured transmission for each 
bolometer may lead to large uncertainties in the estimates of the 
bandpass mismatch leakage (Planck Collaboration VIII 2016), 
so it is difficult to make an accurate prediction of the leakage for 
compact sources. Instead we use the FFP8 simulations (Planck 
Collaboration XII 2016) to assess the effect of the leakage. 

Two sets of FFP8 maps have been generated. One set was 
simulated with the measured bandpasses for each HFI bolome¬ 


ter, so it contains the bandpass mismatch leakage. The other set 
was generated using the average frequency channel bandpass for 
all bolometers in a channel. In this idealized case, there is no 
mismatch between the bandpasses, so no leakage is produced. 
We compare the polarized flux densities of sources from the two 
sets of maps. A suitable quantity for assessing the size of the 
leakage is the difference between the polarized flux density from 
the maps with the leakage, ^bpm, and that from the maps with¬ 
out the leakage, ^nobpm, normalized by the uncertainty due to the 
noise, cr, 

^ _ (‘^bpm “ 6'nobpm) 

cr 

Figure C.l shows this quantity plotted against the input polar¬ 
ized flux density of the source for the 100-353 GHz channels. 
The size of the effect depends on the differences between the in¬ 
dividual bolometer bandpasses and the average frequency chan¬ 
nel bandpass; they are smallest for the 217 GHz channel. For all 
channels the effect of the leakage is smaller than 1 cr for most 
sources. The mean value of A gives the average bias on the po¬ 
larized flux density measurements. It is smaller than 0.06 cr for 
all four channels. Therefore we conclude that the effect is small 
and can be safely ignored. 


C.2.1. Evaluation of HFI polarization measurements 

We use the Q and U maps that have been corrected for the 
leakage of the diffuse temperature components into polarization. 
This does not correct for any leakage from temperature due to 
the compact sources themselves. However, as shown above, this 
effect is subdominant when compared with the statistical uncer¬ 
tainty of the measurements. 

The polarized flux density estimator, P, is evaluated using 
Eq. (1) and the maximum-likelihood estimates, Q and U, are 
extracted from the corresponding Q and U maps by fixing the 
position to that found from the intensity map and strictly assum¬ 
ing that the source is point-like. This selection of the estimator is 
justified, because maximum-likelihood estimators produce min¬ 
imum variance and unbiased estimates, which follow a Gaussian 
distribution with a variance given by 

^ (C.3) 

where N~^ is the inverse power-spectrum of the background 
of the patch-map, and is the beam transfer function as 
function of the bidimensional spatial frequency vector, rj. The 
maximum-likelihood estimator employed for the extraction of 
the Q and U signals was Powell Snakes (PwS; Carvalho et al. 
2012, 2009). The PwS likelihood assumes the background is a 
realization of a homogeneous Gaussian random held, where a 
power-spectrum is known, which is a good assumption for small, 
flat patches cut from the Q and U maps. The PwS package has 
been extensively tested and used inside the Planck Collaboration 
and is known to deliver robust and accurate estimates (Planck 
Collaboration VII2011; Planck Collaboration VIII2011; Planck 
Collaboration XXIX 2014). The flux densities in the Q and U 
maps can be either positive or negative. In order to reduce sys¬ 
tematic effects, which could be induced by the ancillary steps 
of the likelihood evaluation (such as the estimation of the back¬ 
ground power-spectrum of the patch-map), we perform the same 
estimation procedure on both the patch-map and its negative, 
and the average of the two flux-density estimates is taken. This 
procedure helps in stabilizing the polarization signal, especially 
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Fig. C.2. Three different cases of polarized flux density poste¬ 
rior distributions (Eq. C.5). The red curve shows a non-detection 
{p < 1.21), the blue curve shows a marginal detection (1.21 < 
p < 3.0), and the green curve shows a significant detection 
{p > 3.0). 


when tackling regions with complex backgrounds. Equivalent 
estimates of Q and U, and their uncertainties, can be obtained us¬ 
ing an aperture-photometry estimator. The aperture-photometry 
estimator is robust to deviations from the likelihood data model - 
e.g. extended sources, background deviations from Gaussianity, 
or variations in the beam shape - at the cost of slightly larger 
error bars. Eor both estimators the criterion for the acceptance 
of a putative detection was that P had to be 99.99 % signifi¬ 
cant with regard to the null-hypothesis, which is well described 
by a Rayleigh distribution when the assumption is made that 
(Tq o-jj, and that these errors are uncorrelated. In the case 
of acceptance, the polarization angle estimate, is evaluated 
using Eq. (2). The Q and U uncertainties are propagated onto 
the P and ij/ estimates. This is done by assuming the uncertain¬ 
ties are normally distributed and the error bars are small com¬ 
pared with the measured quantities. This approximation holds 
very well given the high significance threshold we have chosen. 
In the case where a putative detection is rejected because it does 
not reach the required significance, we provide the 99 % upper 
limit. 


C.3. Methods to determine polarization properties of 
marginal detections 

Eor the four polarization-sensitive HEI channels, we provide an¬ 
other set of polarized flux-density and polarization angle esti¬ 
mates for sources with marginal detections of polarization. The 
extraction of the polarization signal from the Q and U maps fol¬ 
lows the same procedure described above and in Sect. 2.1. We 
then proceed by assessing the probability of obtaining a given 
measurement of polarization, given a true value of polariza¬ 
tion P. 

Assuming ctq ^ cru ~ cr, the probability of drawing p = 
^Jq^ -f u^jcr given the true value of the polarized flux density 


Pq = Polcr, is 

£ip I cr) = Pr(p I po, o-) = p exp lo(ppo), (C.4) 

where Io(x) is the zero-order modified Bessel function of the 

first kind (Herranz et al. 2012). Using Bayes’ theorem, the pos¬ 
terior distribution of po (the true polarized flux density) given 
the measurement p is 

Pr(Po I P) exp lo(ppo) ^(Po), (C.5) 

where 0(/?o) is the Heaviside step function. We have used the 

Heaviside step function as a prior, since the flux-density param¬ 
eter is intrinsically positive, and by taking the asymptotic form 
of Iq(ppo) it may be shown that po acts as a location parameter 
for large p. Examples of three different posterior distributions 
for the polarized flux density are shown in Eig. C.2. 

The polarization angle distribution was derived using Eq. (2). 
It has been assumed that U and Q are independent Gaussian- 
distributed random variables with means pu = Pq = where 
u, q are maximum-likelihood estimates of U and 2, and cru and 
(Tq are given by Eq. (C.3). Then using the equality for changing 
variables in a probability distribution Pr(^) dO = Pr(f) df, where 
^ = f//2, the probability distribution for 0 is 

Pr)^) = -2 cos(26')“ V (- tan(26')), (C.6) 

where f(w) is a function defined in Hinkley (1969, equation 1). 
The best-fit value is the mode of the distribution and the asym¬ 
metric error-bars were computed using the 95 % highest proba¬ 
bility density (HPD) region of the posterior distribution. HPDs 
are discussed by Box & Tiao (1992) as a general method for 
compressing the information contained within a probability dis¬ 
tribution. Each HPD is uniquely defined by the amount of prob¬ 
ability it encloses and is constructed such that there exists no 
probability density value outside the HPD that is greater than 
any value contained within it. In other words the 95 % HPD con¬ 
tains 95 % of the total probability under the posterior distribu¬ 
tion, such that there is no point outside of this area with a higher 
probability than any point inside the area. This approach allows 
us to provide a best-fit value and the asymmetric 2 cr error-bars 
for the marginal polarization entries in the catalogue. 

We shall refer to this set of measurements as the 
PowellSnakes marginal polarization (PwSPOL) data set. 
PwSPOL permits a proper statistical characterization of fainter 
polarization signals and is therefore able to provide a deeper and 
more complete catalogue without any loss of reliability, as is 
shown in Sect. 3.5. An additional benefit of PwSPOL is that 
it provides a qualitative assessment about an even fainter pop¬ 
ulation of polarized sources which could be valuable as tar¬ 
gets for follow-up observations. It does this by splitting the 
non-detections into two separate groups: clear non-detections, 
where the polarized flux density posterior peaks at no-signal, 
and marginal detections, where the posterior does not peak at no¬ 
signal, but this possibility is still inside the 95 % HPD region of 
the posterior distribution. In Pig. C.3 we compare the histograms 
of the polarized flux-density measurements in the PCCS2 found 
with each approach. The distribution of polarized flux densities 
in the 353 GHz channel does not follow the same pattern as the 
other HEI channels. This is due to the much shallower complete¬ 
ness (see Sect. 3.5) of this channel as compared with the oth¬ 
ers. Although reliable, the polarized flux-density measurements 
of this marginal set are expected to be biased high as result of 
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Fig. C.3. Normalized histograms of the polarized flux density 
from the PCCS2 catalogues for 100 to 353 GHz. The blue line 
shows the detections obtained using the common procedure used 
by both LFI and HFI, whereas the marginal detections, which 
are produced only for HFI, are shown in green. The number of 
sources in each group is shown in the top right of each panel. 


Eddington-type bias and are therefore not suitable for any statis¬ 
tical analysis. 
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